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ABSTRACT 


The  results  of  an  experimental  investigation  on  the  interaction  of  vortex  rings  with  a 
free  surface  are  presented.  Laminar  vortex  rings  with  Reynolds  numbers  in  the  range  2x10^ 
to  10x10^  were  formed  underwater  at  several  incidence  angles  to  the  free  surface.  When  the 
vortex  core  closest  to  the  surface  interacts  with  it,  vortex  lines  are  found  to  disconnect  and 
reconnect  with  the  free  surface  resulting  in  open  vortex  lines  beginning  and  terminating  at 
the  surface.  Under  certain  conditions  a  second  reconnection  process  is  observed  resulting 
in  a  vortex  line  topology  consisting  of  two  semicircles  propagating  away  from  each  other 
on  divergent  paths.  This  double  reconnection  of  voticx  rings  at  a  free  surface  is 
documented  in  detail  by  flow  visualization  and  surface  velocity  field  measurements.  A  new 
technique  was  developed  to  measure  the  velocity  field  at  the  free  surface  using  Laser 
Speckle  Photography  and  automated  digital  image  analysis.  Test  results  show  the  accuracy 
of  the  velocity  measurement  to  better  than  1%. 

It  is  further  shown  that  the  inter^tion  geometry  when  the  vortex  ring  moves 
initially  parallel  to  the  surface  is  independent  of  the  ring  Reynolds  number.  The  effect  of 
Reynolds  number  on  reconnection  time  is  also  investigated  and  showed  to  be  independent 
of  Reynolds  number.  It  is  however  quite  sensitive  to  the  incidence  angle. 

The  effect  of  free  surface  contamination  is  investigated  at  an  incidence  angle  of  90*. 
The  results  show  the  formation  of  secondary  and  tertiary  rings  with  a  contaminated  surface. 
With  a  cleaner  surface  the  formation  of  a  secondary  ring  is  delayed. 
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CHAPTER  I 


INTRODUCTION 

There  is  considerable  interest  on  the  interaction  of  tiirbulent  shear  flow  with  a  free 
surface.  An  important  aspect  of  these  problems  is  the  characterization  of  the  interaction  of 
the  free  surface  with  underwater  vordcity  associated  with  the  large  scale  vortical  structures 
in  the  turbulent  flow.  A  flow  visualization  study  of  the  interaction  of  a  turbulent 
axisymmetric  jet  with  a  free  surface  by  Bernal  and  Madnia  (1988,  1989)  showed  a  complex 
pattern  of  surface  waves  and  surface  nwdons  caused  by  the  vortical  modons  in  the  jet  flow. 
Distinct  dark  spots  were  observed  in  shadowgraph  flow  visualization  images  of  the 
interacdon  region.  These  dark  spots  are  associated  with  vonex  lines  terminadng  at  the  free 
surface  (Berry  and  Hajnal,  1983,  Sterling  et  al.  1987).  Sarpkaya  and  Henderson  (1985) 
studied  the  surface  signature  of  trailing  vordces  generated  by  a  delta  wing.  They  observed 
the  formadon  of  surface  scars  (features  organized  along  the  direcdon  of  the  wing  modon) 
and  striauons  (features  with  general  orientadon  perpendicular  to  the  wing  modon)  on  the 
surface  associated  with  the  vortex  motion.  Dark  spots  indicadng  vortex  lines  terminating  at 
the  free  surface  were  observed  along  the  scars.  Also  when  the  trailing  vortices  underwent 
Crow  instability  (Crow, 1970)  they  observed  the  frxmadon  of  crude  vtxtex  rings  with 
vortex  lines  terminadng  at  the  free  surface.  An  important  consideration  in  the  underwater  jet 
and  underwater  delta  wing  wake  problems  is  that  all  the  voiticity  in  the  flow  is  introduced 
underwater.  As  the  flow  evolves  vordcity  finds  its  way  to  the  free  surface.  The  flow 
process  responsible  for  the  advecdon  of  vordcity  to  the  free  surface  and  generation  of 
vordcity  norma)  to  the  surface  are  not  understood.  The  present  study  of  the  interaaion  of  an 
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axisymmetric  laminar  vortex  ring  with  the  free  surface  was  undertaken  to  gain  some  insight 
into  these  flow  processes.  The  vortex  ring  is  one  of  the  simplest  flow  conflgurations  in 
which  these  processes  can  be  investigated  in  a  fully  three-dimensional  flow. 

A  schematic  diagram  of  the  flow  geometry  with  deflnition  of  parameters  is  shown 
in  Figure  1.1.  The  vortex  ring  is  characterized  by  the  circulation,  F,  the  ring  diameter,  a, 
and  the  core  size,  5.  The  effect  of  vorticity  distribution  within  the  ring  on  the  speed  of 
propagation  has  been  discussed  by  Lamb  (1945),  Saffman  (1970)  and  Widnall  and 
Sullivan  (1973)  among  others.  These  parameters  arc  determined  by  the  formation  process. 
The  effect  of  various  parameters  in  the  formation  process  on  the  circulation,  diameter, 
speed  of  propagation  and  stability  of  vortex  ring  has  been  investigated  by  Saffman  (1975), 
Maxworthy  (1972,1977),  Widnall  and  Sullivan  (1973),  Didden  (1977,1979)  and  Glezer 
(1988)  over  a  wide  range  of  the  Reynolds  number.  It  seems  reasonable  to  expect  that  these 
vonex  ring  parameters  will  influence  the  interaction  with  the  free  surface. 

Cerra  and  Smith  (1983)  repon  a  limited  number  of  observations  in  the  case  of  a 
vortex  ring  impacting  normal  to  the  free  surface.  They  compare  these  results  with  the 
results  of  similar  experiments  in  which  rhe  free  surface  is  replaced  by  a  solid  surface.  They 
found  that  the  evolution  of  the  flow  in  both  cases  is  very  similar.  This  evolution  is 
characterized  by  generation  of  opposite  sign  vorticity  at  the  surface  which  subsequently 
rolls-up  into  secondary  and  even  tertiary  vortex  rings.  The  mechanisms  of  vorticity 
generation  at  the  surface  is  not  identified  although  they  attribute  this  behavior  to  surface 
tension  effects.  It  should  be  noted  that  in  Cerra  and  Smith  experiments  the  Froude  number 
was  small.  Dahm  et  al.(1989)  investigated  the  interaction  of  a  vortex  ring  with  a  density 
interface  at  normal  incidence.  They  identify  the  proper  scaling  parameters  and  show  that  in 
the  limit  of  strong  interface  or  weak  vortex  the  deformation  of  the  interface  is  small  yet 
substantial  opposite  sign  vorticity  is  generated  by  baroclinic  effects.  They  also  show 
striking  similarity  of  this  limit  case  with  the  solid  wall  case.  The  interaction  of  vonex  pairs 


with  the  free  surface  at  normal  incidence  was  investigated  by  Barker  and  Crow  (1977)  and 
more  recently  by  WiUmarth  et  al  (1989)  at  high  Froude  numbers. 

There  is  no  information  on  the  nature  of  the  interaction  of  a  vortex  ring  with  a  free 
surface  at  incidence  angles  other  than  normal  incidence.  Sarpkaya  and  Henderson  (1985) 
proposed  that  vortex  lines  terminating  at  the  surface  are  generated  by  a  process  similar  to 
the  vortex  linking  observed  in  the  Crow  instability  (Sarpkaya,  1983).  The  basis  for  this 
argument  is  that  in  the  case  of  small  Froude  number,  the  free  surface  can  be  approximated 
by  a  flat  nondeformable  surface  orequivalendy  a  symmetry  plane.  Thus,  in  this  limit  the 
flow  field  can  be  expected  to  be  similar  to  the  flow  encountered  when  two  vortex  lines 
interact  across  a  symmetry  plane.  For  vortex  rings  this  type  of  interaction  has  been 
investigated  by  several  authors  including:  Kambe  and  Takao  (1971),  Fohl  and  Turner 
(1975),  Oshima  and  Asaka  (1977),  Schatzle  (1987)  and  Oshima  and  Izutsu  (1988). 

One  striking  phenomenon  encountered  in  these  interactions  is  the  vortex 
reconnection  process  found  when  the  vortex  rings  come  in  contact  along  intersecting  paths. 
Vonex  lines  in  the  incoming  rings  forming  two  separate  closed  loops  disconnect  and 
reconnect  to  form  a  single  closed  loop.  Under  certain  conditions,  as  the  reconnected  vortex 
lines  move  apart,  two  parts  of  the  distorted  vortex  lines  come  into  contact  again  resulting  in 
a  second  reconnection  process  and  the  formation  of  two  vortex  rings  moving  along 
diverging  paths  on  a  plane  perpendicular  to  the  initial  incidence  plane.  This  interaction 
process  is  referred  to  as  a  double  reconnection.  The  double  reconnection  process  was  first 
reported  by  Kambe  and  Takao  (1971)  who  used  smoke  flow  visualization  of  vortex  rings 
in  air  propagating  parallel  to  each  other.  Double  reconnection  was  found  to  occur  for  a 
narrow  range  of  vortex  ring  circulations.  Fohl  and  Turner  (1975)  studied  the  collision  of 
two  identical  rings  directed  toward  one  another  along  intersecting  axes  in  water.  They 
found  a  critical  incidence  angle  for  the  double  reconnection  process  to  occur.  This  critical 
incidence  angle  was  attributed  to  a  resonance  effect  and  they  showed  experimental  results 
for  one  value  of  the  Reynolds  number.  A  systematic  study  of  the  effect  of  the  initial  speed 
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of  vortex  rings  and  distance  between  rings  on  the  interaction  of  rings  moving  side  by  side 
in  air  was  conducted  by  Oshima  and  Asaka  (1977).  Depending  on  the  initial  speed,  they 
classified  the  interactions  into  three  typical  patterns,  single  reconnection,  double 
reconnection  and  multiple  reconnection.  Oshima  and  Izutsu  (1988),  and  Schatzle  (1987) 
have  measured  the  velocity  field  during  the  interaction  of  colliding  vortex  rings.  Oshima 
and  Izutsu  report  experimental  results  of  the  time  dependent  vorticity  field  on  two 
orthogonal  planes  using  x-type  hot  wires  for  two  vortex  rings  moving  side  by  side  in  air. 
They  conclude  that  each  filament  in  the  vortex  tubes  individually  reconnects  and  moves 
away  quickly  in  the  direction  normal  to  the  plane  of  the  filament  due  to  its  strong  curvature. 
Schatzle  studies  the  reconnection  process  by  measuring  the  velocity  field  on  the  symmetry 
plane  with  a  laser-Doppler  velocimeter  fw  two  vortex  rings  colliding  at  a  shallow  angle  in 
water.  He  argues  the  reconnection  process  results  from  vorticity  cancellation  due  to  the 
production  of  out-of-plane  strain. 

It  is  apparent  that  velocity  field  information  is  necessaiy  to  properly  understand 
these  complex  vortical  interactions.  Laser  speckle  velocimetry  (LSV)  or  laser  speckle 
photography  (LSP)  is  a  non-invasive  velocity  measurement  technique  capable  of  measuring 
the  flow  field  in  a  planar  section  at  an  instant  Stetson  (1975),  Vest  (1979)  and  Lautebom 
and  Vogel  (1984)  described  this  technique  in  detail.  Also  Dudderar  et  al  (1986)  published 
an  excellent  review  of  the  development  of  laser  speckle  velocimetry.  Barker  and  Foumey 
(1977),  Grousson  and  Mallick  (1977)  and  Dudderar  and  Simpkins  (1977)  published  almost 
simultaneously  successful  demonstrations  of  LSV  for  Poiseuille  flow.  From  that  time, 
many  researchers  have  been  interested  in  the  LSV  technique,  but  manual  determination  of 
the  velocity  required  lots  of  time  and  lacked  accuracy.  Meynart  (1982,1983)  developed  the 
digital  image  technique  to  analyze  the  data.  Adrian  (1984)  and  Adrian  and  Yao  (1985) 
distinguished  between  the  terms  'speckle  velocimetry'  and  'particle  image  velocimcffy' 
depending  on  the  seeding  density  and  discussed  the  optimum  tracer  particle  size.  An 
important  limitation  of  this  technique  is  that  it  has  a  180 '  ambiguity  in  the  direction  of  the 
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measured  velocity  vector.  Lourenco  et  al  (1986)  and  Adrian  (1986)  have  used  image 
shifting  techniques  to  resolve  the  directional  ambiguity  and  to  increase  the  velocity  dynamic 
range.  The  role  of  photographic  parameters  in  speckle  technique  are  reported  by  Lourenco 
and  Krothapalli  (1987). 

In  the  present  study,  the  interaction  of  a  vortex  ring  with  a  free  surface  is 
investigated.  Experiments  were  conducted  at  several  vortex  ring  generation  conditions  and 
incidence  angles.  The  effect  of  surface  contamination  was  also  explored.  The  conditions 
resulting  in  the  double  reconnection  process  were  determined  and  the  flow  field  of  the 
surface  was  measured  using  LSP.  In  Chapter  2.  the  experimental  apparatus  and 
instrumentation  are  described.  The  details  of  laser  speckle  photography  as  implemented  in 
this  investigation  and  the  automatic  digital  image  processing  technique  used  are  discussed 
in  (Chapter  3.  The  results  of  the  experin^nts  for  the  various  incidence  angles  and  the  detail 
flow  field  measurements  are  presented  and  discussed  in  Chapter  4.  A  summary  of  the  main 
results  and  conclusions  of  the  investigation  is  given  in  Chapter  S. 


CHAPTER  II 


FLOW  FACILITY  AND  INSTRUMENTATION 


2J _ The.raalUt 

The  experiments  were  conducted  in  the  flow  facility  shown  schematically  in  Hgure 
2.1.  The  facility  consists  of  a  free  surface  water  tank  and  a  vortex  generator  assembly.  The 
free  surface  water  tank  is  made  of  glass  76.2cm  wide,76.2cm  high  and  167.6cm  long.  The 
tank  was  filled  with  Ann  Arbor  city  water  filtered  through  a  5  |im  sediment  removal  filter. 

Preliminary  experiments  showed  large  sensitivity  of  the  interaction  process  to  the 
presence  of  surface  acdve  agents  in  the  water  as  discussed  in  Section  4.2.  Several  facility 
elements  were  designed  and  incorporated  to  minimize  the  amount  of  surface  active  agents  in 
the  water.  As  shown  in  Figure  2.1,  the  entire  facility  was  covered  with  plastic  sheet  to 
prevent  dust  panicles  and  other  contaminants  fnxn  settling  on  the  water  surface.  Also  a 
stand-up  pipe  connected  to  the  drain  was  used  to  control  the  water  level  in  the  tank.  Fine 
adjustment  of  the  vortex  generator  depth  was  obtained  by  small  adjustments  of  the  length  of 
the  pipe.  To  insure  a  clean  surface  a  continuous  flow  of  water  was  allowed  into  the  tank 
which  was  removed  through  the  stand-up  pipe.  This  technique  preferentially  removes  water 
from  the  free  surface  retixiving  with  it  the  surface  active  material  before  it  can  accumulate. 

The  condnuous  flow  of  water  was  set  up  at  least  one  day  before  the  experiments 
were  conducted  and  throughout  the  duradon  of  the  experiments.  This  flow  was  interrupted 
only  a  few  minutes  before  the  generadon  of  each  vortex  ring  to  allow  any  current  in  the 
tank  to  die  out  Under  these  condidons  the  water  temperature  in  the  tank  was  maintained  at 
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8.5*  C  ±  0.2*  C  for  the  duration  of  the  experiments.  The  surface  tension  of  water  in  the 
tank  was  measured  using  the  capillary  tube  technique.  Typical  values  measured  when  these 
procedures  were  followed  were  in  the  range  71-73  dynes/cm  which  arc  within  5%  of  the 
value  for  pure  water  at  that  temperature. 

L2 _ Vortex  Generator 

A  photograph  of  the  vortex  generates  used  in  the  investigation  is  shown  in  Figure 
2.2.  As  shown  in  Figure  2.3,  the  vortex  generator  consisted  of  a  inlet  flow  diverter  (0.3  cm 
thick  and  5  cm  in  diameter),  a  foam  section  (2.5  cm  thick),  a  honeycomb  section  (3.8  cm 
long  and  0.4  cm  in  diameter)  and  a  contraction  section.  The  contraction  section  consisted  of 
a  circular-arc  shaped  nozzle  (17.8  cm  radius)  with  a  contraaion  ratio  of  16  to  1.  The  exit 
orifice  had  circular  cross-section  2.54  cm  in  diameter. 

The  vortex  generator  was  mounted  on  a  support  bar  which  in  turn  was  mounted  on 
a  unismit  structure  affixed  to  the  top  of  the  tank,  as  shown  in  Figure  2.4.  The  vortex 
generator  coarse  depth  adjustment  was  provided  by  the  depth  adjustment  bar.  The  vortex 
ring  incidence  angle,  a,  was  adjusted  by  rotating  the  vortex  generator  around  a  pivot  point 
at  the  end  of  the  depth  adjustment  bar.  A  threaded  hole  and  lock  nut  were  used  to  secure  the 
vortex  generator  at  the  desired  angle. 

The  flow  system  used  to  generate  the  vonex  rings  is  also  shown  in  Figure  2.4.  This 
system  consists  of  a  pressurized  water  reservoir,  a  solenoid  valve  and  associated  control 
components.  In  a  typical  operation  water  in  the  reservoir  tank  is  pressurized  with  shop  air 
using  the  pressure  regulator.  The  vtxtex  ring  is  generated  by  opening  the  solenoid  valve  for 
a  shon  time.  The  parameters  that  control  die  circulation  strength  of  the  resulting  vortex  ring 
are  the  reservoir  pressure,  the  water  needle  valve  setting  and  the  solenoid  valve  (^)en  time 
duration.  A  submersible  solenoid  valve  (ASCO  8031B43, 4  psi  pressure,  1.9  cm  orifice 
size)  was  mounted  just  behind  the  vtxtex  generator  (see  Figure  2.2)  to  reduce  oscillation 


after  the  solenoid  valve  was  closed.  A  stainless-steel  reservoir,  10  cm  in  diameter  by  42  cm 
long,  was  connected  to  the  solenoid  valve  through  the  needle  valve  (1.27  cm  orifice)  and  a 
flexible  tube  for  setting  the  arbitrary  angle  of  the  vortex  generator.  The  pressurized  air  was 
supplied  to  the  reservoir  through  the  pressure  regulator  and  needle  valve  which  were 
mounted  on  the  tq)  of  the  reservoir.  The  reservoir  and  vortex  generator  were  filled  with 
water  supplied  through  the  top  of  the  reservoir,  or  by  the  suction  of  the  air  inside  the 
reservoir  through  the  shut-off  valve.  To  eliminate  air  bubbles  inside  the  vortex  generator, 
water  was  warmed  up  to  room  temperature  and  used  to  fill  the  vortex  generator  in  a  vertical 
position.  The  temperature  of  the  water  inside  the  generator  was  allowed  to  equilibrate  with 
the  water  temperature  in  the  free  surface  tank  by  running  tap  water  continuously  in  the  tank 
through  the  night  with  a  rubber  plug  at  the  exit  orifice  of  the  vortex  generator. 

The  solenoid  valve  opening  time  was  controlled  by  the  timing  circuit  shown  in 
Figure  2.5.  It  consists  of  a  timer  chip  and  a  solid-state  relay  to  energize  the  valve. 
Additional  timing  components  were  used  as  shown  in  Figure  2.5  to  synchronize  the 
generation  of  the  vortex  ring  with  other  equipment  used  for  flow  visualization  as  described 
below. 


_ Flow  Visualization 


Several  flow  visualizadon  techniques  were  used  in  this  investigation.  Simultaneous 
flow  visualizadon  of  the  free  surface  and  of  the  cross-section  of  the  vonex  ring  were 
obtained  using  the  system  shown  schemadcally  in  Figure  2.6.  The  surface  deformadon 
induced  by  the  vortex  ring  was  visualized  using  a  shadowgraph  technique.  A  collimated 
beam  was  formed  by  posidoning  a  mercury  arc  lamp  at  the  focal  point  of  a  300  mm 
diameter  f/5  spherical  mirror.  The  collimated  beam  was  directed  perpendicular  to  the 
undisnirbcd  free  surface  by  means  of  a  30  cm  x  30  cm  first  surface  flat  minor.  The 
refracted  pattern  caused  by  deformadon  of  the  free  surface  was  viewed  on  a  screen 
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positioned  at  an  angle  of  45*  with  respect  to  the  free  surface.  Simultaneous  visualization  of 
the  cross-section  of  the  vortex  rings  was  obtained  by  adding  a  fluorescent  dye  (Rhodamine 
6G)  at  a  concentration  of  3  ppm  to  the  water  in  the  vortex  generator.  The  flow  was 
illuminated  with  a  laser  light  sheet  The  sheet  was  formed  from  the  green  (514  nm)  line  of 
an  Argon-Ion  laser  by  means  of  a  collimator  and  a  cylindrical  lens  (15  cm  in  focal  length). 
The  sheet  was  positioned  along  the  symmetry  plane  of  the  ring  perpendicular  to  the  free 
surface.  The  thickness  of  the  sheet  in  the  region  of  interest  was  approximately  1  mm.  Both 
the  vortex  ring  cross  section  and  surface  deformation  were  simultaneously  recorded  on  the 
same  photographic  film  or  video  image. 

Experiments  using  the  simultaneous  visualization  of  the  cross-section  view  and  free 
surface  deformation  (See  Section  4.2. 1.1)  showed  that  both  images  could  not  be 
simultaneously  focused  on  the  film.  While  slightly  out  of  focus  images  could  be  tolerated  in 
the  relatively  lower  resolution  video  images,  it  was  not  acceptable  for  higher  resolution 
photographic  work.  A  better  quality  image  of  the  flee  surface  deformation  was  obtained 
separately  using  the  imaged  shadowgraph  system  shown  schematically  in  Figure  2.7.  In 
addition  to  the  set  up  for  shadowgraph  technique,  another  first  surface  flat  mirror  was 
installed  on  top  of  the  flee  surface  tank  to  direct  the  beam  into  another  300  mm  in  diameter 
f/5  spherical  mirror  and  a  lens  located  at  the  focal  point  of  the  spherical  mirror.  This  lens 
focused  the  shadow  image  formed  43  cm  above  the  flee  surface  on  the  film  plane.  Different 
magnifications  of  the  image  were  obtained  using  different  focal  length  lenses. 

The  photographs  of  the  vonex  ring  cross  section  were  obtained  using  a  Nikon  35 
mm  camera  equipped  with  a  motor  drive  controlled  by  the  timing  circuit  shown  in  Figure 
2.5.  The  solenoid  valve  and  motor  drive  timers  were  triggered  at  the  same  time.  The  motor 
drive  adjustable  timer  setting  determines  the  time  delay  between  the  formation  of  the  vonex 
ring  and  the  instant  at  which  the  picture  was  obtained.  A  135  mm  fl2.8  lens  was  used  with 
a  31  mm  extension  ring  to  provide  proper  magnification.  The  camera  was  located  close  to 
the  tank  to  capture  the  whole  evolution  of  vonex  ring  during  the  interaction  with  the  free 
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surface.  The  photographs  of  the  free  surface  deformation  using  the  imaged  shadowgraph 
system  were  obtained  using  the  same  camera  body  with  a  lens  as  indicated  above. 

2J — Hot  Film  Velocity  Measurement 

Hot  film  velocity  measurements  were  conducted  to  determine  the  circulation 
strength  of  the  vortex  ring  and  the  formation  characterisdcs  of  the  vortices.  A  constant 
temperature  hot  film  anemometer  with  a  standard  TSI  quartz  coated  cylindrical  senscn*  was 
used.  The  sensca*  fiequency  response  is  specified  by  the  manufacturer  as  80  KHz.  The 
sensor  axis  was  positioned  perpendicular  to  the  exit  orifice  along  the  centerline.  A  block 
diagram  of  the  electronics  and  data  acquisition  system  is  shown  in  Figure  2.8.  A  constant 
temperature  anemometer  circuit  (TSI  model  1750)  was  connected  to  the  decade  resistance 
box  (TSI  model  1056)  with  a  resolution  of  0.01  Q.  The  output  of  the  anemometer  was  DC 
shifted  and  amplified  using  a  Tektronix  AM  501  operational  amplifier.  To  digitize  the 
output  of  the  difierential  amplifier,  a  LeCroy  Model  8210  Transient  Digitizer  was  used  at  a 
sampling  rate  of  IKHz.  The  hot  film  probe  was  mounted  on  a  computer  controlled  traverse 
mechanism  (Willmarth,  1977).  A  Transicoil  Model  U-217094  Gearhead  Stepper  Motor 
was  used  to  move  the  probe. 

The  probe  calibration  was  performed  using  the  vratex  generatOT  and  the  same  free 
surface  tank  facility.  The  probe  was  located  at  the  vortex  generate^’  exit  and  the  exit  velocity 
was  determiMd  by  measuring  the  time  required  to  discharge  a  fixed  volume  out  of  the 
reservoir.  The  vortex  generator  exit  velocity  was  varied  by  changing  the  reservoir  pressure. 
For  each  vortex  generator  exit  velocity,  the  output  voltage  from  the  amplifier  was  recorded. 
A  calibration  curve  was  obtained  by  a  fourth  order  polynomial  fit  to  the  measured  data.  The 
estimated  error  of  the  calibration  curve  is  1.7%  of  the  maximum  exit  velocity.  To  minimize 
drift  of  the  velocity  signal,  the  water  temperature  was  maintained  within  0.2*  C  during  die 
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calibration  and  the  experiment.  Also  the  water  was  passed  through  a  softener  filter  in 
addition  to  the  5  }xm  sediment  removal  filter. 


CHAPTER  III 

FREE  SURFACE  VELOCITY  FIELD  MEASUREMENT 

3J^_  General  Consideration 

In  these  experiments  the  free  surface  velocity  field  was  measured  using  Laser 
Speckle  Photography.  The  fundamental  concept  is  to  recrad  on  the  same  photogr^hic  plate 
two  time-displaced  images  of  the  flow.  The  flow,  in  this  case  the  free  surface,  is  seeded 
with  small  tracer  particles.  The  fluid  velocity  component  on  a  plane  parallel  to  the 
photographic  plate  is  determined  from  the  displacement  of  the  tracer  panicles  during  the 
time  between  exposures.  The  implementation  of  the  technique  can  be  divided  into  three 
separate  steps:  (1)  Photographic  Recording,  (2)  Flow  Image  Analysis  and  (3)  Velocity 
Field  Analysis.  While  these  three  steps  are  interrelated  from  the  point  of  view  of  accuracy, 
resolution  and  overall  efficiency  of  the  technique,  they  represent  three  physically  di^erent 
procedures  using  different  equipment  which  are  better  discussed  separately.  After  a  brief 
description  below,  step  (2)  is  discussed  in  detail  in  Section  3.2,  step  (3)  is  discussed  in 
Section  3.3  and  step  (1)  is  discussed  in  Section  3.4. 

In  step  (1)  a  photographic  image  of  the  flow  seeded  with  particles  is  obtained.  The 
details  of  the  technique  used  for  the  free  surface  velocity  measurement  are  described  in 
Section  3.4.  Test  images  were  also  obtained  and  used  for  development  and  calibration  of 
the  software  used  in  steps  (2)  and  (3).  For  this  purpose  specklegrams  were  obtained  using 
the  rotating  disk  set  up  shown  in  Figure  3.1.  A  lucite  disk  mounted  on  a  dc  motor  was 
rotated  at  an  angular  velocity  of  1000  rpm.  A  pulsed  ruby  laser  with  a  typical  exposure  time 
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of  20  ns  was  used  to  illuminate  the  disk.  Double  exposures  were  obtained  with  a  dme 
between  exposure  of  60  ps.  The  resulting  speckle  patterns  were  recorded  on  a  holographic 
plate,  AGFA-Gevaert  10E75,  which  has  a  resolving  power  of  better  than  3000  lines/mm. 

In  step  (2)  the  displacement  of  the  particles  is  measured  by  directing  a  narrow 
collimated  beam  of  coherent  light,  typically  an  unexpanded  He-Ne  laser  beam,  through  the 
photographic  plate  at  the  point  of  interest  Young's  fringes  are  formed  by  interference  of 
the  diffraction  patterns  of  the  displaced  speckles  and/cH*  particle  images  on  the  film.  These 
interference  fringes  are  aligned  perpendicular  to  the  direction  of  the  surface  displacement 
The  fringe  image  formed  at  the  back  focal  plane  of  a  lens  were  analyzed  using  a  digital 
image  analysis  system.  The  displacement  of  the  panicles  between  exposures,  d,  is  related 
to  the  fringe  spacing,  df,  by  the  equation 


(3.1) 


where  X  is  the  wavelength  of  the  laser  used  to  form  the  fringes,  fp  is  the  focal  length  of  the 


transforming  lens  and  M  is  the  magnificadon  of  the  imaging  system  used  to  take  the 
photograph.  Then  the  local  flow  velocity  is  given  by  the  reladon 


(3.2) 


where  At  is  the  dme  between  exposures. 

A  typical  fringe  image  is  shown  in  Figure  3.2.  Because  the  fringe  image  is  formed 
by  coherent  light  a  random  speckle  pattern  is  superposed  on  the  fringes.  The  characterisdc 
scale  of  this  speckle  pattern  is  determined  by  the  interrogadon  beam  diameter.  It  is  apparent 
that  the  Young's  fringes  will  not  form  if  the  particle  displacement  on  the  photographic  plate 
is  larger  than  the  interrogation  beam  diameter.  In  the  present  system  this  limitation  implies 
that  as  the  flow  displacement  increases  the  fringe  spacing  decreases  until  it  becomes  of  the 
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order  of  the  speckle  size  and  can  no  longer  be  measured.  The  technique  is  also  limited  by 
the  fact  that  the  fnnge  spacing  can  not  be  measured  when  it  is  comparable  to  the  size  of  the 
diffraction  halo  around  the  intenogadon  beam.  This  limit  corresponds  to  the  displacement 
of  the  particles  on  the  film  being  small  compared  to  the  resolution  of  the  imaging  system 
since  it  is  these  small  features  on  the  photographic  plate  that  determine  the  size  of  the 
diffraction  halo.  An  important  limitation  of  the  technique  is  that  there  is  a  180*  ambiguity  in 
the  direction  of  the  measured  velocity  vector.  A  velocity  bias  technique  can  resolve  this 
ambiguity  (Lourenco  et  al  (1986)  and  Adrian  (1986))  and  also  improve  the  dynamic  range 
of  the  technique.  In  the  present  measurements  velocity  bias  was  not  used.  The  direction  of 
the  velocity  vector  was  determined  from  global  analysis  of  the  flow  field.  The  velocity  field 
was  measured  over  the  region  of  interest  by  successively  scanning  the  photograph  and 
measuring  the  fringe  spacing  and  orientation  on  a  point  by  point  basis. 

Finally  in  step  (3),  the  measured  velocity  vector  field  is  analyzed  to  determine 
relevant  pmper^es  like  the  component  of  the  vonicity  vector  normal  to  the  plane  of  the 
measurement  and  the  deformation  rate  tensor  components  on  the  plane  of  the  measurement. 

_ Fringe  Image  Analysis 

3.2.1  Image  Processing  System 

The  fringe  spacing  and  direction  on  the  fringe  image  were  measured  using  a  high 
speed  Digital  Image  Processing  System.  A  schematic  diagram  of  the  Image  Processing 
System  is  shown  in  Figure  3.3.  A  CW  He-Ne  laser  (Jodon  Inc.  Model  #  HN-IOG)  is  used 
to  interrogate  the  flow  image  at  a  point  The  fringe  image  is  formed  in  the  back  focal  plane 
of  a  Fourier  transform  lens.  A  beam  stop  was  positioned  at  the  center  of  the  fringe  image  to 
block  the  bright  spot  due  to  the  undiffracted  light  A  host  computer  (2[EOS  286)  controls 
the  image  processor  and  a  two  component  traverse  system  (Newport  Corp.  Model  8SSc 
Programmable  Controller  and  Model  SSO-series  Actuators).  A  video  camera  (DAGE-MTI 
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Inc.  Model  65)  is  used  to  digitize  the  fringe  image  into  a  512x512  pixel  anay  with  8  bits 
density  resolution  for  a  total  of  1/4  Megabyte.  The  resulting  digital  image  is  stored  and 
analyzed  in  a  Gould  FD5000  Image  Processing  System.  The  FD5000  contains  two 
channels  of  image  memory.  Each  channel  is  divided  into  four  512x512  pixel  banks.  Each 
pixel  is  stored  in  memory  as  a  16  bit  value.  The  upper  8  bits  of  each  pixel  are  used  for 
graphic  overlays  or  additional  image  memory.  The  lower  8  bits  give  the  image  density  in 
the  range  from  0,  dark,  to  255,  white.  A  real  time  Feedback  Processor  (FBP)  is  the  main 
component  in  the  FD5000.  The  layout  of  the  FBP  is  shown  in  Figure  3.4.  The  input  to  the 
FBP  is  selected  ftom  one  of  the  image  memories  or  a  predefined  constant.  The  input  data 
can  be  modiHed  using  a  lookup  table  (LUT).  The  LUTs  are  loaded  from  the  computer  and 
replace  each  input  value  with  a  corresponding  value  in  the  table.  Data  is  combined  in  the 
arithmetic  logic  unit  (ALU)  which  performs  simple  arithmetic  and  logic  functions.  Data  is 
then  sent  to  an  image  memory  directly  ot  by  way  of  an  output  LUT.  The  accumulator  can 
be  used  to  sum  the  intensities  of  pixels  over  the  entire  image.  The  image  input  to  the  FBP 
can  be  modified  by  the  zoom/scroll  device.  This  device  controls  a  window  on  the  image 
memory  and  can  move  the  window  over  the  image  memory.  In  addition,  the  FD5(X)0  has  a 
cursor  device  and  trackball  which  allow  the  user  to  select  points  on  the  image. 

3.2.2  Image  Enhancement 

As  shown  by  Figure  3.2,  the  fringes  on  the  fringe  image  appear  superposed  on  a 
speckle  pattern.  From  the  point  of  view  of  measuring  the  fringe  spacing  this  speckle  pattern 
is  noise  and  must  be  removed.  An  image  analysis  technique  was  developed  to  filter  the 
speckle  noise.  It  is  convenient  to  describe  this  technique  on  the  basis  of  the  spectral  content 
of  the  fringe  image.  The  power  spectral  density  (PSD)  along  the  horizontal  (x  axis)  and 
vertical  axis  (y  axis)  were  obtained  using  a  modified  Welch  method  (Welch,  1967).  For  the 
PSD  along  the  horizontal  direction,  the  whole  image  was  first  sectioned  into  512  segments 
of  512  pixels  length  along  the  horizontal  direction.  Each  segment  corresponds  to  a 
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horizontal  line  of  this  image.  A  modified  periodogram  for  each  segment  was  computed 
using  the  Fast  Fourier  Transform  algorithm.  Then  the  sp>ectral  estimate  was  calculated  from 
the  average  of  these  512  periodograms.  The  PSD  along  the  vertical  direction  was  also 
obtained  using  the  same  method.  In  this  case  each  segment  corresponds  to  a  vertical  line  of 
the  image.  The  measured  horizontal  and  vertical  PSDs  for  the  fringe  image  in  Rgure  3.2 
are  shown  in  Figure  3.5.  Throughout  this  discussion  the  PSD  curves  are  normalized  for 
unit  area  under  the  curve.  The  PSD  curves  show  the  main  features  of  the  fringe  image  that 
must  be  considered  in  the  analysis.  Rrst,  there  is  considerable  energy  at  low  wavenumbers 
which  corresponds  to  background  nonuniformity  on  the  original  image.  Second,  on  each 
PSD  there  is  a  distinct  spike  which  corresponds  to  the  fringe  pattern.  Third,  there  is 
significant  broadband  noise  associated  with  the  random  speckle  pattern.  It  is  the  objective 
of  the  image  enhancement  technique  to  isolate  the  second  feature,  the  fringe  pattern,  from 
the  other  two  which  are  measurement  noise. 

To  eliminate  the  random  speckle  noise  while  retaining  the  fringe  pattern,  the  image 
was  averaged  along  the  expected  fringe  direction  characterized  by  a  direction  angle  Os.  The 
direction  of  the  averaging  was  set  manually  using  two  cursors  on  the  TV  monitor  and  the 
trackball.  Figure  3.6  shows  the  scheme  for  line  averaging.  The  image  was  shifted  along  the 
desired  direction  using  the  scroll  device  in  the  FD5(XX)  and  stored  in  another  channel. 
Corresponding  pixels  on  the  original  image  and  the  shifted  image  were  added  and  the  result 
divided  by  two  using  the  FBP.  The  process  was  repeated  N  times.  For  each  successive 
iteration  the  magnitude  of  the  image  shift  was  increased  by  a  factor  of  2.  Using  this  scheme 
an  average  length  of  2^  pixels  was  implemented  with  only  order  N  operations.  A 
photograph  of  the  resulting  line  averaged  image  for  N=5  is  shown  in  Figure  3.7.  The 
corresponding  PSDs  arc  shown  in  Figure  3.8.  Both  the  line  averaged  image  and  the  PSD 
results  show  that  the  speckle  noise  was  successfully  removed  while  the  fringes  were 
retained.  This  image  is  still  contaminated  with  considerable  background  nonuniformity. 
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To  eliminate  the  background  nonunifonnity  the  original  fringe  image  was  low-pass 
filtered  and  the  resulting  background  image  was  subtracted  from  the  line  averaged  image. 
The  low-pass  filter  was  implemented  using  an  equal  weight  square  window  of  2^  x  2^ 
pixels.  The  same  averaging  scheme  used  for  the  line  average  was  utilized  which  results  in 
2N  operations.  A  value  of  N=5  was  found  adequate  to  generate  a  good  background  image. 
A  photograph  of  the  image  after  background  subtraction  is  shown  in  Figure  3.9  and  the 
results  for  the  corresponding  PSDs  are  shown  in  Figure  3. 10.  This  image  and  the  PSD 
show  that  the  background  nonuniformity  was  removed  The  spikes  in  the  PSD  associated 
with  the  fiinges  are  more  pronounced  in  this  image  with  a  signal  to  noise  ratio  of  the  order 
of  100.  The  width  of  these  spikes  however  is  relatively  large.  To  further  sharpen  the  spikes 
contrast  enhancement  was  applied  to  this  image.  Intensity  levels  below  120  were  set  to  0 
and  above  140  were  set  to  255.  Intensity  levels  between  120  and  140  were  redistributed 
linearly  from  0  to  255.  This  process  increased  the  visibility  of  the  fringes  by  sharpening  the 
intensity  transitions  associated  with  the  fringes.  A  photograph  of  the  contrast  enhanced 
image  is  shown  in  Figure  3.1 1.  The  corresponding  PSDs  are  shown  in  Figure  3.12.  This 
image  is  now  ready  for  analysis  by  the  autocorrelation  technique  described  below.  The 
width  of  the  spikes  in  the  PSD  is  of  the  order  of  the  resolution  (1/512  cycles/pixel)  of  the 
measurements.  This  implies  that  the  spatial  frequency  of  the  fringes  can  only  be  determined 
within  that  uncertainty.  This  corresponds  to  an  uncertainty  of  the  fringe  spacing 
measurement  of  the  order  of  20%  for  the  case  shown  in  Figure  3. 1 2.  The  autocorrelation 
technique  described  below  was  developed  to  reduce  this  uncenainty  to  less  than  1%. 

3.2.3  Autocorrelation  Technique 

The  fringe  spacing  and  orientation  were  determined  from  the  autocorrelation 
function  of  the  intensity  distribution  on  the  contrast  enhanced  image  along  two 
perpendicular  axes.  The  normalized  autocorrelation  along  each  axis  was  computed  from  the 
relation 
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where  L  =  200  pixels.  Typical  autocorreladon  results  along  two  perpendicular  axes  are 
shown  in  Figure  3.13.  The  first  zero  crossing  in  this  figure  occurs  at  one  fourth  of  the 
fringe  spacing.  Subsequent  zero  crossings  are  spaced  by  one  half  of  the  fringe  spacing.  It 
should  be  noted  that  the  amplitude  of  the  fluctuation  of  the  correlation  function  is  large  for 
offsets  of  the  order  of  200  pixels  (almost  half  image).  This  is  a  direct  consequence  of 
improved  visibility  of  the  fringes  in  the  contrast  enhanced  image  (Figure  3.1 1). 

The  orientation  of  the  fringes,  af,  was  determined  by  measurement  of  the  fringe 

spacing  along  two  perpendicular  directions.  The  accuracy  of  the  orientation  measurement  is 
improved  by  selecting  the  two  axes  as  follows.  If  the  expected  fringe  angle,  Og,  was  within 
-  3/8  Ji  <  <  -  jc/8  or  ji/8  <  aj  <  3ic/8,  then  axes  parallel  to  the  x  axis  and  the  y  axis 

were  used.  If  as  was  in  the  region  as  -  jc/2  <as  <  -  3x/8  or  -  x/8  <  ttj  <  7t/8  or  37t/8  <as 
<  7t/2,  then  axes  oriented  at  45*  and  -  45*  with  respect  to  x  axis  were  chosen.  This  choice 
insures  that  the  fringe  direction  is  not  parallel  to  one  of  the  measurement  directions. 

The  fringe  spacing  along  each  direction  was  determined  from  the  first  and  last  zero- 


crossings  of  the  autocorrelation  function  within  a  200  pixel  window.  These  zero-crossings 
were  determined  directly  using  an  efficient  algorithm  (number  of  operations  « log2 
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(offset)).  From  the  location  of  the  first  zero  crossing  Zi  and  final  zero  crossing  zf,  the  fringe 
spacing  was  determined  as 


ks  Int 

where  k  is  the  number  of  half  cycles  of  the  autocorrelation  function  between  zj  and  zf.  The 
value  of  k  was  verified  since  it  must  be  even  for  a  negative  slope  of  the  autocorrelation 
function  at  Zf  and  odd  for  a  positive  slope  of  the  autocorrelation  function  at  zf.  The 
uncertainty  on  the  df  mcasunement  is  determined  by  the  uncertainty  on  zf  which  is  ±  0.5 
pixels.  Since  Z(  is  of  the  order  of  2(X)  pixels.  The  relative  error  in  df  is  less  than  1  %. 

3.2.4  Sensitivity  Study 

An  important  possible  limitation  of  this  technique  is  that  it  requires  an  input 
parameter  namely  the  expected  fringe  direction  0$.  A  sensitiyity  study  was  conducted  to 
determine  the  yariation  of  measured  fringe  spacing  and  measured  fringe  direction  as  a 
function  of  the  expected  fringe  direcdon.  The  results  for  the  fringe  spacing  and  fringe 
direction  are  shown  in  Figure  3.14  and  Figure  3.15  respectiyely  for  a  typical  image.  The 
expected  fringe  angle  was  yaried  oyer  a  range  of  ±  5*.  The  measured  fringe  spacing  dm  is 
within  1%  of  the  actual  fnnge  spacing  df.  The  measured  fringe  angle  ocm  is  within  ±  0.4*  (± 
0.01  rad)  of  the  fringe  angle  af.  The  results  suggest  that  if  the  change  of  fnnge  direction 
between  two  consecutiye  points  measured  on  the  same  flow  image  is  small  (within  5*),  the 
yalue  of  the  fringe  direction  measured  at  the  preyious  point  can  be  used  as  the  expected 
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fringe  angle  for  the  present  point  In  this  case,  the  expected  fringe  direction  angle  was 
needed  only  for  the  first  point  The  entire  flow  image  can  be  analyzed  without  operator 
interaction. 


3.2.5  Implementation 

A  computer  program  was  developed  to  implement  the  flow  field  velocity 
measurement  using  this  technique.  A  flow  chart  for  this  program  is  given  in  Rgure  3. 16. 
The  program  runs  in  the  host  computer,  see  Rgure  3.3  and  interacts  with  the  Image 
Processing  System.  This  program  automatically  scans  the  flow  image  with  the  traverse 
system.  It  can  operate  in  automatic  mode  and  manual  mode.  In  automatic  mode  there  is  no 
operator  interaction  after  the  iiutial  set  up.  In  manual  mode  the  expected  fringe  direction  at 
each  point  must  be  entered  by  the  operator.  To  help  this  process  the  fringe  direction  of  the 
previous  point  is  displayed  on  the  current  point  fringe  image.  The  results  of  the 
measurements  along  with  the  position  information  are  stored  on  disk  flies  for  subsequent 
plotting  and  analysis. 

The  specklegram  obtained  in  the  rotating  disk  experiment  was  analyzed  with  this 
program  operated  in  automatic  mode.  The  velocity  vector  was  measured  in  a  square  matrix 
of  32  X  32  points  spaced  1  mm,  for  a  total  of  1024  measurement  points.  The  elapsed  time 
required  for  the  measurement  was  6  seconds/poinL  A  significant  amount  of  this  time  (2.5 
seconds4>oint)  was  spent  traversing  the  specklegram.  The  elapsed  time  per  point  is  larger  in 
manual  mode  because  of  the  time  required  to  input  the  expected  fnnge  angle. 

_ Velocity  Field  Analysis 

Several  flow  properties  including  the  vorticity  field  and  several  components  of  the 
deformation  rate  tensor  can  be  calculated  from  these  data.  Before  calculating  these  other 


properties  a  vector  plot  of  the  velocity  field  was  obtained.  Figure  3.17  is  the  vector  plot  of 
the  velocity  field  in  the  rotating  disk  experiment  The  data  was  measured  in  a  square  region 
on  the  disk  31  mm  x  31  mm.  The  accuracy  of  data  throughout  this  region  is  better  than  1%. 
Because  the  disk  is  in  solid  body  rotation  these  data  can  also  be  used  to  estimate  the  errors 
associated  with  the  determination  of  the  vorticity  and  deformation  rate  fields. 

The  component  of  the  vorticity  vector  normal  to  the  plane  of  the  measurement,  O), 
was  determined  using  the  relation 
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where  A  is  the  area  of  a  surface  element  around  the  point  of  interest  and  the  integral 
represents  the  circulation  around  that  surface  element  The  value  of  co  at  each  point  was 
calculated  by  considering  a  surface  element  bounded  by  the  8  neighboring  points.  The  line 
integral  along  those  points  was  evaluated  using  Simpson's  rule.  For  the  rotating  disk 
experiment  the  disk  is  in  solid  body  rotation  so  that  the  vorticity  at  a  point  equals  twice  the 
angular  velocity  and  is  uniform  across  the  region.  From  the  data  shown  in  Figure  3. 17  the 
vorticity  was  evaluated  throughout  the  region.  Because  the  vorticity  is  uniform  in  this  case 
a  statistical  measure  of  the  accuracy  of  the  vorticity  measurement  is  obtained  from  the  mean 
and  standard  deviation  of  the  values  in  the  sample.  The  mean  vonicity  for  the  sample  is  o)  = 
210.67  sec**  and  the  standard  deviation  is  48.9  sec**.  The  mean  value  is  in  good  agreement 
with  the  expected  value  of  209.44  sec**  derived  from  the  angular  velocity  of  the  disk  1000 
rpra  The  scatter  around  this  value  is  consistent  with  the  expected  uncertainty  of  the 
measurement  given  by  the  ratio  of  the  velocity  measurement  uncertainty  divided  by  the 
distance  between  measured  points  (i.e. »  50  sec**  in  this  case). 

From  the  velocity  field  measurements  the  components  of  the  deformation  rate 
tensor,  £xx,  £yy  and  £xy ,  on  the  plane  of  the  measurement  can  be  calculated.  These 


c(Hnponents  of  the  deformation  rate  tensor  were  calculated  using  a  second-<xder  central 
difference  scheme.  Thus  at  a  point  in  the  sample  (xi,  yO  the  spatial  derivatives  were 
calculated  using  the  expression 
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where  /ij  can  be  either  the  u  or  v  component  of  the  velocity  at  the  corresponding  point. 

The  deformation  or  strain  rate  information  is  better  expressed  in  terms  of  the 
directions  of  the  principal  axes  and  the  principal  strain  rate  values  along  those  axes.  The 
Jacobi  method  was  used  to  determine  one  principal  axis  direction  with  respect  to  the  x  axis 
9  as 
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with  the  positive  sign  used  if  £xx  >  £yy  and  the  negative  sign  otherwise.  The  magnitude 
of  the  strain  rate  along  6  is  given  by 


€ll  *  Cxx  cos^O  +  £xy  sin20  +  £yy  sin^G. 


(3.11) 


Of  course  the  other  principal  axis  is  perpendicular  to  the  9  direction  and  the  corresponding 
magnitude  of  the  strain  rate  is  given  by 


£22  “  £xx  sin29  -  £xy  sin29  +  £yy  cos29. 


(3.12) 
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In  the  case  of  the  rotating  disk  data  shown  in  Figure  3. 17,  where  the  disk  is  in 
solid  body  rotation,  all  components  of  the  deformation  rate  tensor  are  zero  everywhere  on 
the  disk.  Again  the  mean  and  standard  deviation  of  the  computed  values  of  £xx,  £yy  and 

£xy  give  a  statistical  measure  of  uncertainty  for  these  parameters.  The  mean  and  standard 
deviation  for  £xx,  £yy  and  £xy  were  evaluated  for  the  rotating  disk  sample.  Typical 
values  for  the  mean,  £xx,  and  standard  deviation,  £xx' ,  of  the  normal  deformation  rate 
along  the  x  axis  are  £xx  =  -8.56  sec**,  £xx’  =  97.88  sec**  respectively. 

These  results  are  consistent  with  the  expected  uncertainty  of  the  measurement 

2A _ Photographic  Recording 

The  laser  speckle  photography  technique  was  used  to  measure  the  velocity  field  at  a 
free  surface.  Preliminary  measurements  were  conducted  on  the  region  of  interaction  of  a 
tubulent  underwater  jet  with  the  free  surface.  Typical  results  will  be  discussed  in  this 
section.  The  technique  was  later  applied  to  study  the  interaction  of  a  vortex  ring  with  the 
free  surface.  A  schematic  diagram  of  the  facility  and  optical  system  is  shown  in  Figure 
3. 18.  For  the  jet  experiment  a  small  jet  tank  (described  in  Appendix  A)  located  inside  the 
glass  tank  was  used  to  produce  the  jet  flow.  The  jet  centerline  was  located  one  diameter 
(2.54  cm)  below  the  free  surface  and  the  jet  exit  velocity  was  25  cm/sec.  Tracer  panicles 
which  followed  the  fluid  motion  at  the  surface  were  added  to  the  free  surface  tank.  A 
Copper  Vapor  Laser  was  used  to  illuminate  the  particles.  The  optical  system  was  similar  to 
the  one  used  fix  the  surface  imaging  shadowgraph  flow  visualization.  A  53  cm  focal  length 
lens  was  used  which  resulted  in  a  magnification  of  0.4.  The  jet  exit  velocity  used  was  kept 
low  to  minimize  the  free  surface  deformation  thus  reducing  the  refraction  effects  on  the 
image. 

The  surface  tracers  used  in  the  experiments  were  glass  microballoons.  The  density 
of  the  microballoons  was  0.15  g/cc  and  the  mean  diameter  was  approximately  45  pm.  The 


concentration  of  tnicroballoons  at  the  free  surface  was  5x10^  cm*^.  A  photograph  of 
microballoons  is  shown  in  Figure  3.19.  The  size  variation  of  the  tnicroballoons  in  this 
picture  is  typical  of  the  samples  used  in  the  experiments.  Because  of  the  small  size  of  the 
microballoons  they  are  expected  to  cause  very  small  deformations  of  the  surface.  An 
important  consideration  is  the  ability  of  tracers  to  stay  at  the  free  surface  and  follow  the 
surface  modon.  The  Stokes  characterisdc  time  of  the  microballoons  in  water  is  8  ^s  so  they 
follow  the  velocity  fluctuadons  well.  The  limitadon  of  the  technique  as  implemented  in  this 
experiment  is  the  fact  that  the  jet  fluid  was  not  seeded  with  tracer  particles.  Hence  regions 
in  the  free  surface  where  the  fluid  originated  from  the  jet  may  not  have  been  properly 
seeded  with  tracer  pardcles  and  velocity  measurement  may  not  be  possible  there. 

In  order  to  properly  expose  the  film  a  Copper  Vapor  Laser  was  programmed  to 
produce  two  bursts  of  light  each  consisdng  of  five  pulses.  The  time  between  pulses  in  each 
burst  was  100  |is.  The  film  used  to  obtain  the  flow  images  was  high  resolution 
photographic  plate  AGFA-Gevaert  10E56.  The  dme  between  bursts  was  5  ms.  It  is  during 
this  dme  that  the  significant  displacement  of  the  tracer  pardcles  occur.  The  value  was 
determined  from  consideration  of  the  expected  flow  velocity,  the  speckle  size  and  the 
intertogadon  beam  diameter. 

A  typical  fringe  image  is  shown  in  Figure  3.20.  In  addition  to  the  Young's  fringes, 
elliptical  fringes  are  also  observed.  Because  the  spacing  between  these  ellipdcal  fringes  is 
comparable  to  the  Young's  fringe  spacing,  the  digital  image  processing  program  could  not 
be  used  in  some  regions.  A  vector  plot  of  the  surface  flow  velocity  field  is  shown  in  Figure 
3.21.  The  frame  area  in  this  plot  corresponds  to  2<xJDo<9  and  -2<y/Do<2,  where  y  is  the 
coordinate  axis  perpendicular  to  the  flow  direcdon  on  the  plane  of  the  free  surface.  At  the 
jet  centerline  the  surface  flow  velocity  increases  with  downstream  distance  as  a  result  of  the 
interaction  of  the  jet  with  the  free  surface.  The  maximum  velocity  measured  in  this  region  is 
15.6  cm/s.  In  contrast  the  velocity  away  from  the  centerline  is  quite  low  and  the  minimum 
velocity  measured  here  is  0.7  cm/s.  The  entrained  flow  moves  in  the  upstream  direcdon  as 


it  approaches  the  vortical  region.  Also  the  velocity  vector  rotates  as  the  fluid  enters  the 
vortical  region.  The  component  of  the  vordcity  vector  perpendicular  to  the  free  surface  was 
evaluated  from  the  measured  velocity  field  using  the  same  method  applied  to  the  rotating 
disk  experiment.  A  contour  plot  of  the  vorticity  is  shown  in  Figure  3.22.  The  uncertainty  of 
measurement  of  (i)D(AJe  is  estimated  as  ±  0.07.  As  expected  the  entrained  flow  is  by  and 
large  irrotational.  Vorticity  at  the  free  surface  is  introduced  by  the  interaction  of  the  surface 
with  the  vortical  flow  under  it.  The  scale  of  the  vortical  features  on  the  free  surface  is  of  the 
order  of  the  jet  width.  These  data  clearly  show  the  importance  of  the  large  scale  structure  of 
the  Jet  in  the  interaction  with  the  free  surface. 

To  investigate  the  interaction  of  the  vortex  ring  with  the  free  surface,  the  speckle 
photographs  were  obtained  with  the  Copper  Vapor  Laser  in  the  optical  configuration  shown 
in  Figure  3.23.  The  vortex  generator  was  inclined  20  ‘  with  respect  to  the  free  surface.  The 
depth  of  the  center  of  the  vortex  gun  exit  cross  section  was  h/Do  =  1.72  below  the  free 
surface.  AGFA-Gevaert  10E56  photographic  plate  was  positioned  parallel  to  the  free 
surface  1  cm  above  the  free  surface  to  record  the  specklegram.  The  Copper  Vapor  Laser 
was  programmed  to  produce  two  bursts  of  light  each  consisting  of  5  pulses  synchronized 
with  the  timing  circuit  that  operated  the  solenoid  valve  (Figure  2.5).  The  time  between  burst 
was  set  at  6  ms.  The  tracer  particles  used  for  this  experiment  were  the  same  microballoons 
as  used  for  jet  interaction  experiment  The  surface  contaminants  at  the  free  surface  were 
removed  by  a  continuous  surface  current  as  described  above.  Prior  to  the  experiment  the 
surface  current  was  turned  off.  After  the  water  came  to  rest  tracer  particles  were  deposited 
uniformly  on  the  free  surface  by  tapping  a  sieve  which  contained  the  particles  above  the 
region  of  interest  When  tracer  particles  stopped  moving,  the  vortex  ring  was  generated  and 
the  speckle  photograph  was  taken  after  the  proper  time  delay.  When  specklegrams  were 
processed  using  the  digital  image  processing  program,  different  focal  length  lenses  (53.3 
cm  and  30.5  cm)  were  used  depending  on  the  displacement  of  the  speckle  patterns  on  the 
specklegram  to  improve  the  velocity  dynamic  range.  The  results  are  discussed  in  Chapter  4. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 

_ The  Characteristics  of  the  Vortex  Rings 

Experiments  were  conducted  to  determine  the  characteristics  of  vortex  rings 
obtained  with  the  vortex  generator.  In  these  experiments  the  vonex  rings  were  formed  in  an 
effectively  unbounded  fluid  by  locadng  the  generator  a  long  distance  away  from  the  free 
surface  and  the  walls  of  the  facility.  The  main  parameters  of  interest  are  the  propagation 
speed,  Ut,  the  circulation  strength,  F,  and  the  vorticity  distribution  within  the  vortex.  Of 
these  paranwters  the  propagation  speed  and  circulation  are  relatively  easy  to  measure.  The 
vorticity  distribution  on  the  other  hand  is  very  difficult  to  determine.  It  is  convenient  to 
introduce  two  parameters  associated  with  the  vorticity  distribution  in  the  ring;  the  vortex 
ring  core  diameter,  a,  and  the  core  size,  6,  as  the  simplest  characterization  which  still 
retains  the  most  important  dynamics  of  the  ring. 

Hve  experimental  conditions  are  considered  in  this  investigation.  They  were 
selected  based  on  the  data  by  Glezer  (1981,1988)  and  after  an  extensive  flow  visualization 
study.  In  all  these  cases  the  vortex  rings  were  laminar  at  formation  and  remained  laminar 
over  the  region  of  interest.  The  various  conditions  are:  (1)  the  reservoir  pressure  p=2  psig, 
needle  valve  opening  N=3/4  turns  and  solenoid  valve  opening  time  To=0.2  sec.  (2)  p=l.S 
psig,  N=3/4  turns,  To=0.6  sec.  (3)  p=4  psig,  N=3/4  turns,  To=0.2  scc.(4)  p=  4psig,  N=1 
turn,  To*0.2  sec.  Finally,  (5)  p=4  psig,  N=3/4  turns,  To=0.6  sec. 
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To  obtain  the  propagation  speed  and  circuladon  of  the  vortex  ring,  hot  film  velocity 
measurements  were  conducted  at  several  locations  along  the  axis  of  symmetry  of  the 
generator.  A  typical  velocity  time  history  for  condidon  (4)  is  shown  in  Figure  4.1.  The  hot 
film  probe  was  located  5  nozzle  exit  diameters  from  the  generator  exit  plane.  The  time 
origin  corresponds  to  the  opening  of  the  solenoid  valve.  It  can  be  seen  that  the  velocity  is 
very  small  up  to  1.7  seconds  after  the  generation  of  the  vortex  ring.  The  flow  velocity  is 
maximum  at  2.15  seconds  after  generation  and  then  decreases.  The  flow  velocity  is 
effectively  zero  after  2.7  seconds.  Similar  velocity  traces  were  obtained  at  several 
downstream  locations  and  for  the  selected  vortex  ring  generation  conditions.  For  the 
purpose  of  determining  the  velocity  of  propagation  of  the  vonex  ring,  the  arrival  time  of 
ring  at  the  probe  location  is  defined  as  the  time  at  which  the  velocity  is  maximum.  A  plot 
of  the  position  of  the  vortex  ring  as  a  function  of  arrival  time  is  presented  in  Figure  4.2. 
The  data  at  each  operating  condition  is  well  approximated  by  a  straight  li  .e.  A  least  square 
fit  to  the  data  was  used  to  determine  the  speed  of  propagation  of  the  vortices  given  in  Table 
4.1.  The  uncertainty  of  the  measurement  is±  3%.  The  propagation  speed  was  also 
measured  from  video  recordings  of  the  evolution  of  the  vortices.  The  results  give  a  slightly 
higher  value  than  the  one  obt'^ined  by  hot  film  measurements.  This  discrepancy  is  attributed 
to  the  uncertainty  of  the  video  system  frame  rate. 

The  circulation  strength  was  determined  from  the  velocity  measured  along  the 
centerline  using  two  different  techniques.  First,  the  velocity  measured  at  the  center  of  the 
exit  orifice  on  the  exit  plane  was  used.  In  order  to  determine  the  circulation  it  is  necessary 
to  assume  a  uniform  velocity  profile  at  the  exit  except  for  a  thin  boundary  layer  at  the  wall. 
The  circulation.  To.  is  determined  from  the  total  flux  of  vorticity  through  the  exit  plane  as 
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(4.1) 


where  y  is  the  radial  coordinate  measured  from  the  wall,  Do/2  is  the  value  of  y  at  the  orifice 
center,  Uo(t)  is  the  velocity  at  the  centeriine  and  Tq  is  the  duration  of  the  water  pulse  (Note: 
it  is  assumed  U(t)  =  0  for  t  <  0  and  t  >  Tq).  The  values  of  To  determined  from  the  velocity 
traces  measured  at  the  exit  plane  of  the  vortex  generator  are  presented  in  Table  4. 1 . 

The  second  technique  used  to  determine  the  circulation  is  based  on  the  velocity 
traces  measured  downstream  along  the  axis  of  the  flow.  At  an  instant  the  circulation  of  the 
ring  is  given  by. 


'  =  ^  u  •  dl , 


(4.2) 


where  the  contour  includes  the  ring  axis  of  symmetry  and  is  closed  at  infinity.  If  the 
circulation  can  be  assumed  constant  as  the  vortex  ring  convects  pass  a  fixed  downstream 
location  then 


and  therefore 


dl  =  U,dt, 
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r=u,ru(t)dt. 
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(4.4) 


The  values  of  F  calculated  using  this  expression  are  ploned  in  Figure  4.3  as  a  function  of 
downstream  position.  These  results  show  that  the  measured  circulation  of  the  ring 
decreases  with  downstream  distance  for  cases  (2)  and  (4)  while  is  constant  for  the  other 
two  cases.  In  all  cases  T  reaches  a  constant  value  further  downstream.  The  decrease  of  the 
measured  circulation  with  downstream  distance  indicates  that  the  vonex  ring  generation 
process  is  no  completed.  Downstream  of  X/Dq  =  3  the  measured  circulation  is  constant 


which  indicates  a  well  formed  vortex.  The  measured  values  of  F/v  after  fonnation  are  given 
for  all  conditions  in  Table  4.1.  The  uncertainty  of  these  values  is  estimated  as  4%. 

Equation  4.1  is  based  on  the  simple  model  of  the  generation  process  which 
characterizes  the  exhaust  of  fluid  flom  the  generator  as  a  cylindrical  slug  of  length  Lq  and 
diameter  Dq  (Saffman,197S).  If  we  define  a  mean  velocity  as 

_  1  fT. 

*  T  *  I  »  r  j.. 


Uo  =  :rf  *Uo(t)dt 
^  O  Jo 


where  Uo(t)  is  assumed  uniform  over  the  orifice  area,  the  equivalent  length  of  the  slug  is 


(t)  dt  =  Uo  T„ 


Following  Glezer  (1981,1988)  the  impulse  associated  with  the  slug  is  given  by 


I  =  pA 
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J  Uo  (t)dt=: 
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where  A  is  the  cross-section  area  at  the  exit  of  the  generator.  The  nondimensional  slug 
length  L^o  is  an  important  parameter  of  the  formation  process  (Glezer,  1988).  The  values 
of  Lo/Dq  were  determined  from  the  velocity  traces  at  the  exit.  The  results  are  given  in  Table 


It  is  apparent  that  the  values  of  the  circulation  obtained  from  the  velocity  measured 
downstream  differs  from  the  values  obtained  for  the  velocity  measured  at  the  exit  plane  of 
the  generator.  Figure  4.4  is  a  plot  of  the  normalized  circulation  of  the  well  formed  vortex 
FA'  as  a  function  of  the  circulation  determined  from  the  generator  exit  velocity  Fq/v. 

Also  plotted  in  this  figure  are  the  results  of  Didden  (1977).  The  results  at  low  values  of 
Fq/v  indicate  an  excess  circulation  of  the  well  formed  vortex  compared  to  the  value 
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determined  firom  the  velocity  at  the  exit  Case  (4)  at  a  larger  value  of  Fq/v  shows  a  lower 
value  of  FA'  compared  to  Fq/v. 

This  discrepancy  between  the  values  of  F  and  Fq  clearly  indicates  that  equation  4.1 
does  not  provide  an  accurate  estimate  of  vortex  ring  circulation.  Further  the  observation  that 
at  some  conditions  Fi>Fo  while  at  other  Fo>F  suggest  that  competing  mechanisms  which 
tend  to  increase  or  decrease  the  circulation  of  the  vortex  ring  are  at  work  during  the 
formation  process. 

In  an  experimental  investigation  Didden  (1979)  showed  that  there  are  two  different 
mechanisms  which  are  not  accounted  for  in  the  slug  model.  Fu^t,  the  velocity  distribution 
across  the  exit  plane  is  not  uniform.  The  velocity  outside  the  internal  wall  boundary  layer  is 
larger  than  at  the  centerline.  This  effect  is  very  strong  during  the  initial  stages  of 
deformation.  This  increased  velocity  is  the  result  of  streamline  curvature  at  the  exit  of  the 
nozzle.  It  follows  that  the  vorticity  flux  is  considerably  larger  than  what  could  be  expected 
from  the  velocity  measured  at  the  centerline.  The  second  mechanism  is  the  flux  of  vonicity 
associated  with  the  boundary  layer  on  the  external  surface  of  the  generator  induced  by  the 
vortex  ring.  This  vorticity  is  of  opposite  sign  compared  to  the  vorticity  in  the  ring  and 
results  in  a  reduction  of  the  circulation  of  the  ring.  It  should  be  noted  that  the  frrst 
mechanism  is  dominant  in  the  initial  stages  of  formation  while  the  second  mechanism  will 
be  stronger  in  the  later  stages  of  formation.  Thus,  it  can  be  expected  that  the  circulation 
ratio  F/Fq  will  depend  on  L^q-  To  examine  the  effect  of  Lo/Do»  the  circulation  ratio  F/Fq 
is  plotted  in  Figure  4.5  as  a  function  of  Lq/Dq  for  the  present  dr  ind  Didden's  (1977) 
data.  It  is  apparent  that  as  the  length  of  the  slug  is  increased  the  circulation  ratio  approaches 
unity.  These  results  are  in  agreement  with  Didden's  data  and  are  consistent  with  Didden's 
(1979)  detailed  study  of  the  formation  process. 

The  initial  vortex  ring  diameter,  a,  was  measured  on  video  recordings  of  only  the 
nozzle  exit  region  using  a  zoom  lens.  The  vortex  core  size,  5,  could  not  be  measured 
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directly.  It  was  computed  from  the  measured  values  of  Ut,  a  and  F  using  the  theoretical 
result  for  a  constant  vorticity  core  (Lamb, 1932) 


The  calculated  values  of  5  are  also  given  in  Table  4. 1 . 


(4.8) 


£2 — Vortex  Ring  Interaction  with  Free  Surface 


The  interaction  of  a  vortex  ring  with  the  free  surface  was  investigated  as  a  function 
of  incidence  angle  a  (See  Figure  1.1).  The  case  of  the  vortex  ring  moving  parallel  to  the 
free  surface  (C£=0*)  was  investigated  first  The  initial  results  suggested  a  significant  impact 
on  the  evolution  of  flow  by  the  presence  of  surface  active  agents.  These  results  as  well  as 
experiments  with  a  cleaner  surface  are  discussed  in  Section  4.2.1.  The  e^ect  of  surface 
contamination  was  investigated  in  more  detail  at  normal  incidence  (a=:90*).  These  results 
are  discussed  in  Section  4.2.2.  The  general  features  of  the  interaction  at  other  values  of  a 
are  discussed  in  Section  4.2.3. 


4.2.1  Vonex  Ring  Parallel  to  the  Free  Surface 

4.2.1. 1  Contaminated  Free  Surface 

Several  experiments  were  conducted  using  the  jet  tank  apparatus  described  in 
Appendix  A  which  generated  vortices  propagating  parallel  to  the  free  surface,  a^O*.  No 
special  care  was  taken  to  reduce  contamination  of  the  free  surface.  All  the  experiments  were 
conducted  with  the  exit  orifice  centerline  located  at  h/Do~  1.2S  below  the  free  surface.  The 
photographs  of  a  typical  evolution  of  the  vortex  ring  at  a  Reynolds  number  Fq/v  =  7.4  x 
103  and  Froude  number  Fo/(ga3)i/2  3  0.41  are  shown  in  Figure  4.6.  In  these  photographs 
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the  cross-sectional  view  and  surface  deformation  are  simultaneously  visualized  (See 
Appendix  A).  Figure  4.6  (a)  shows  the  initial  evolution  of  the  vortex  ring.  A  small 
depression  of  the  free  surface  is  observed  above  the  vortex  ring.  In  Figure  4.6  (b),(c)  the 
upper  part  of  the  vonex  ring  interacts  with  the  free  surface.  As  a  result  of  the  interaction  the 
closed  vortex  lines  in  the  ring  open  and  are  now  terminated  at  the  free  surface.  Figure  4.6 
(d)  shows  that  as  the  lower  part  of  the  vortex  ring  interacts  with  the  free  surface,  opposite 
sign  vorticity  is  generated  at  the  surface  which  rolls  up  in  a  secondary  vonex.  The  vonex 
core  then  rebounds  away  from  the  firee  surface. 

From  video  recordings  of  the  vortex  ring  evolution,  several  measurements  were 
obtained  to  characterize  the  interaction  process.  The  free  surface  influences  the  evolution  of 
the  vonex  ring  causing  a  small  deflection  of  the  vortex  axis  towards  the  free  surface.  The 
extent  of  interaction  was  characterized  by  measurement  of  the  downstream  location  of  the 
upper  vortex  core  when  the  vortex  ring  outline  first  reached  the  free  surfaceXu*  and  the 
downstream  location  of  the  lower  vortex  con  at  the  point  of  closest  proximity  to  the  free 
surface,Li.  The  results  for  Ly/a  and  Li/a  as  a  function  of  a^  for  several  Reynolds  numbers 
are  presented  in  Figure  4.7.  While  these  results  suggest  an  exponential  dependence  of  L/a 
as  a  function  of  a/h,  a  Reynolds  number  effect  can  not  be  ruled  out  from  these  data.  A  more 
detailed  study  of  this  effect  was  conducted  with  a  clean  surface.  Further  discussion  of  this 
result  is  presented  in  Section  4.2. 1.2. 

Detailed  data  on  the  evolution  of  vortex  rings  were  obtained  at  a  Reynolds  number 
of  7.4  X  103  and  Froude  number  0.41.  The  normalized  downstream  location  of  the  vortex 
ring  cores  as  a  function  of  nondimensional  time  are  shown  in  Figure  4.8.  The  path  of  the 
lower  vortex  core  can  be  approximated  by  two  straight  lines  of  different  slope.  The  slope 
change  occurs  at  a  nondimensional  time  of  20.  A  least  square  fit  to  the  data  gives  a 
nondimensional  velocity  Uta/To  =  0.36  before  the  interaction  and  Uta/To  =  0.066  after  the 

interaction.  The  large  reduction  in  propagation  speed  can  be  attributed  to  the  combined 
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effect  of  the  change  in  radius  of  curvature  of  the  vortex  lines  as  well  as  the  clockwise 
vorticity  generated  at  the  free  surface. 

In  Figure  4.9,  the  normalized  distance  of  the  vortex  cores  to  the  free  surface, Y/a, 
and  normalized  distance  between  vortex  cores  measured  on  the  free  surface.  D/a,  are 
plotted.  The  vortex  cores  were  first  observed  on  the  free  surface  at  a  nondimensional  time 
of  13.0,  after  the  vortex  tube  reconnection  process  is  completed.  The  distance  between 
vortex  cores  on  the  free  surface  and  the  distance  of  the  lower  vortex  core  to  the  free  surface 
oscillate  after  the  interaction.  The  two  oscillatory  motions  are  out  of  phase  suggesting  an 
oscillation  of  the  opened  vonex  ring  eccentricity  with  a  period  of  approximately  Totp/a^  = 
50.  This  value  compares  well  with  the  theoretical  calculation  (Widnall  and  Sullivan,  1973) 
for  small  amplitude  oscillations  (mode  0=2)  of  a  thin  cored  vortex  ring  having  the  initial 
circulation  and  the  measured  diameter  and  propagation  speed  after  the  interaction,  Totp/a^  = 
41. 

4.2. 1.2  Clean  Free  Surface 

The  experiments  with  contaminated  surface  suggest  that  vorticity  generation  at  the 
surface  can  play  an  important  role  in  the  interaction  of  a  vortex  ring  with  a  free  surface  at  an 
incidence  angle  of  a=0*.  Similar  experiments  were  conducted  with  the  vortex  generator  and 
a  clean  surface  obtained  as  described  in  Section  2.1.  These  experiments  were  conducted  at 
several  Reynolds  numbers  and  initial  vortex  ring  depths  to  clarify  the  role  of  these 
parameters  on  the  extent  of  the  interaaion. 

The  experiments  with  a  clean  surface  showed  an  evolution  of  the  vortex  rings 
similar  to  the  contaminated  surface  case  (Figure  4.6).  The  interaction  of  the  upper  core  with 
the  free  surface  always  resulted  in  vortex  line  reconnection  as  with  a  contaminated  surface. 
The  main  difference  was  in  the  interaction  of  the  lower  core  with  surface  which  did  not 
result  in  opposite  sign  vorticity  generation.  The  vortex  core  did  not  rebound.  It  appeared  to 
break  down  as  a  result  of  the  interaction.  As  for  the  clean  surface  the  extent  of  the 


interaction  was  determined  from  video  recordings  for  the  flow  evolution  at  three  different 
depths  h/Do  =  1.25, 1.5, 1.75  and  the  five  Reynolds  numbers  (FoA')  in  Table  4.1.  The 
results  are  shown  in  Figure  4.10  using  the  same  notation  as  for  the  contaminated  free 
surface.  These  results  confirm  that  the  normalized  extent  L/a  depends  exponentially  on  a/h, 
and  is  independent  on  Reynolds  number.  Also  comparing  these  results  to  the  results  in 
Figure  4.7  for  a  contaminated  free  surface^t  appears  that  the  interaction  occurs  further 
downstream  for  a  clean  surface  at  low  values  of  a/h. 

4.2.2  Normal  Incidence  to  the  Free  Surface 

The  experiments  at  oc=0*  showed  significant  vorticity  generation  at  a  contaminated 
free  surface.  In  a  related  investigation  Willmanh  and  Hirsa  ( Private  Communication,  1989) 
studied  the  interaction  of  a  vortex  pair  with  the  free  surface,  and  also  found  significant 
effects  attributed  to  the  presence  of  surface  active  agents  at  the  surface. 

To  examine  the  effect  of  surface  active  agent  on  the  interaction  of  vortex  rings  with 
the  surface  the  normal  incidence  (a=:90’)  was  chosen.  This  incidence  angle  was  selected  to 
avoid  vortex  line  reconnection  processes  during  the  interaction,  thus  focussing  these 
experiments  on  vorticity  generation  at  the  free  surface.  The  vortex  generator  was  positioned 
90*  to  the  free  surface  and  the  nozzle  exit  was  at  h/Do=3.5  below  the  free  surface.  For  the 
experiments  with  a  contaminated  free  surface,  the  water  was  left  in  the  tank  overnight 
without  a  surface  cover  and  without  surface  current  flow.  The  clean  surface  was  obtained 
as  described  in  Section  2.1. 

Photographs  of  the  evolution  of  the  vortex  ring  cross-section  were  obtained  at  a 
Reynolds  number  (FqA')  of  3800  and  Froude  number  of  about  0.2  for  both  a  contaminated 
and  a  cleaner  free  surface.  Photographs  are  shown  in  Figure  4.1 1.  Each  photograph  was 
obtained  on  a  different  realization  of  the  flow.  The  location  of  the  free  surface  is  indicated 
on  each  photograph.  Figure  4.1 1  (A)  shows  the  interaction  of  the  vortex  ring  with  a 
contaminated  surface.  During  the  interaction  secondary  and  tertiary  vonex  rings  are  formed 
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having  vorticity  opposite  to  that  of  the  incoming  vortex  ring.  This  evolution  of  the  vonex 
ring  is  very  similar  to  the  evolution  at  a  solid  surface  (Ceira  and  Smith,  1983).  This 
observation  suggests  that  the  water  on  the  free  surface  did  not  move  during  the  interaction. 
Observations  of  small  particles  on  the  free  surface  confirmed  this  conjecture  indicating  that 
capillary  forces  prevented  fluid  motion  at  the  surface.  A  sequence  of  photographs  of  the 
vortex  ring  interaction  with  a  cleaner  free  surface  is  shown  in  Figure  4.1 1  (B).  The 
evolution  of  the  vortex  ring  is  closer  to  what  would  be  expected  in  the  absence  of  capillary 
forces  at  the  surface.  The  self-induced  outward  motion  of  the  vortex  core  persists  for  a 
longer  time.  Also  the  vonex  core  size  is  significantly  reduced  due  to  stretching. 

Figure  4.12  and  4.13  are  plots  of  the  trajectory  of  the  vonex  ring  core  for  the 
contaminated  and  cleaner  surfaces  respectively.  In  the  case  of  a  contaminated  surface,  the 
formation  of  the  secondary  and  tertiary  vonex  rings  limits  the  outward  motion  of  the 
primary  vonex  to  two  and  a  half  times  the  initial  radius.  Also  the  induced  velocity  field  of 
these  vonices  causes  the  primary  vonex  core  to  rebound,  i.e.  it  moves  away  from  the 
surface.  For  the  cleaner  surface  case  shown  in  Figure  4. 13  the  outward  motion  persists  to  a 
radius  of  three  and  a  half  times  the  initial  radius.  Although  there  is  dye  left  behind  the 
primary  vonex  core,  the  amount  of  vorticity  associated  with  this  fluid  should  be  small 
because  the  trajectory  of  the  primary  vonex  core  is  not  altered.  A  secondary  vonex  ring  is 
formed  similar  to  the  one  found  in  the  contaminated  surface  case  which  causes  the 
trajectory  of  the  primary  vortex  to  bounce  away  from  the  surface. 

In  summary  these  results  show  that  surface  contamination  can  have  a  profound 
effect  oh  the  interaction  of  vonex  rings  with  the  surface.  Vorticity  generation  at  the  free 
surface  and  the  subsequent  roll  up  into  secondary  and  tertiary  vonices  causes  a  very  large 
change  of  the  evolution  of  the  ring.  Even  in  the  case  of  very  small  amounts  of  surface 
contaminants,  like  in  the  clean  surface  case,  they  eventually  accumulate  to  generate  a 
secondary  vonex  ring. 


The  observations  at  a  =  0*  showed  vortex  line  reconnection  to  occur  as  a  result  of 
the  interaction  of  the  vortex  ring  with  the  free  surface.  It  was  expected  on  the  basis  of  this 
result  and  existing  evidence  on  the  interaction  of  colliding  vortex  rings  that  the  interaction  at 
an  inclined  incidence  angle  will  also  result  in  vortex  lines  reconnection.  An  extensive  flow 
visualization  study  was  conducted  to  determine  the  general  features  of  the  interaction  in  the 
range  0*  <  a  <  50*  at  vortex  ring  generation  conditions  (1)  to  (4)  in  Table  4.1  and  with  a 
clean  surface.  A  second  objective  of  the  study  is  to  gain  some  understanding  of  the 
reconnection  process.  In  all  cases  considered  the  Froude  number  of  the  interaction, 
Tl{g2?)^ri,  given  in  Table  4.1  is  small  which  indicates  small  surface  deformations  during 
the  interaction. 

Simultaneous  flow  visualization  of  the  underwater  flow  and  of  the  free  surface 
deformation  was  conducted  at  incidence  angles  a  «  5*,  10*,  15*.  20*,  30*,  40*,  45*  and 
50*,  for  several  vonex  generator  depths  and  operating  conditions  (1)  to  (5)  in  Table  4.1. 
The  evolution  of  the  vortices  was  recorded  on  video  tape.  The  results  show  that  the 
interaction  of  the  upper  vortex  core  with  the  surface  always  results  in  vortex  line 
reconnection  up  to  the  maximum  incidence  angle  tested  of  a  =  45*.  The  interaction  of  the 
lower  vortex  core  with  the  surface  was  a  sensitive  function  of  the  vonex  incidence  angle 
and  generation  conditions.  For  incidence  angle  a  ^  5*  the  lower  core  interaction  resulted  in 
vonex  breakdown  at  all  Reynolds  numbers.  For  incidence  angles  in  the  range  40*  £  a  ^ 

45*  there  was  also  vortex  breakdown  after  the  lower  core  interaction  except  the  case  (3) 
(F/v  =  4960).  For  10*  ^  a  ^  30*  cases  (1),(2),  (4)  and  for  10*  ^  a  S  45*  case  (3)  the 
interaction  of  the  lower  core  with  the  free  surface  also  results  in  vonex  line  reconnection.  In 
these  cases  where  the  upper  and  lower  core  reconnect  at  the  surface  the  vonex  ring  splits 
into  two  half  rings  each  having  vonex  lines  starting  and  terminating  at  the  free  surface  (A 
detailed  sequence  of  flow  visualization  pictures  of  this  type  of  interaction  is  described 
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below  in  Section  4.3.1).  These  two  half  vortex  rings  propagate  independent  of  each  other 
after  the  interaction.  The  direction  of  propagation  depended  on  the  flow  conditions.  For 
case  (3)  15*  ^  a  ^  30*  and  case  (4)  15*  S  a  ^  30*  the  half  vortex  rings  propagate  away 
from  each  other  after  the  interaction.  For  case  (1)  15*  <  a  ^  30*,  case  (2)  15*  <  a  <  30* 
and  case  (3)  40*  ^  a  iS  45*  the  half  vortex  rings  propagate  parallel  to  each  other  after  the 
interaction.  For  case  (1)  a  =  10*  the  half  vortex  rings  propagate  towards  each  other  after 
the  formation  is  completed. 

One  important  parameter  of  the  vortex  line  reconnection  phenomena  is  the  time 
required  for  the  process  to  be  completed.  It  could  be  argued  that  the  vortex  reconnection 
time  should  scale  with  the  local  vortex  properties.  On  the  other  hand  global  parameters  like 
the  vortex  ring  incidence  angle,  a,  and  the  Reynolds  number  may  also  play  a  role.  The 
reconnection  time,  tf,  was  defined  for  the  purpose  of  this  study  as  the  time  elapsed  from  the 
initial  contact  of  the  vortex  ring  outline  with  the  free  surface  to  the  time  when  the  vortex 
cores  are  first  observed  as  distinct  dark  spots  on  the  free  surface.  The  measurements  were 
limited  to  the  reconnection  time  of  the  upper  vortex  core.  The  nondimensional  time  Ftr/S^  is 
plotted  as  a  function  of  Reynolds  number  FA'  in  Figure  4.14  for  several  vortex  ring 
conditions.  These  results  show  that  the  nondimensional  reconnection  time  is  independent  of 
the  Reynolds  number.  The  effect  of  incidence  angle  is  significant  and  results  in  an  increase 
of  the  normalized  reconnection  time  with  reduced  incidence  angle.  At  a  =  0*  there  is 
somewhat  more  scatter  in  the  results.  The  effect  of  increasing  initial  depth  is  to  reduce  the 
reconnection  time.  In  any  case  the  reconnection  time  at  a  =  0*  is  longer  than  the 
reconnection  time  for  a  >  0*  in  all  cases  investigated. 

Schatzle  (1987)  argued  that  vorticity  cancellation  during  the  reconnection  can  only 
be  attained  by  viscous  diffusion  and  therefore  the  viscosity  coefficient  must  appear  on  the 
proper  normalization  of  the  reconnection  time.  Our  measurements  do  not  support  this 
conjecture  since  the  vortex  reconnection  time  at  fixed  incidence  angle  is  independent  of 
Reynolds  number.  This  inviscid  scaling  was  proposed  by  Ashurst  and  Meiron  (1987).  The 
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free  surface  velocity  field  measurements  and  detailed  visualizations  described  in  the  next 
section  provide  some  insight  on  this  result 

U - The  Double  Reconnection  Process 

The  double  reconnection  process  was  studied  in  detail  for  vortex  ring  generation 
conditions  corresponding  to  case  (3)  in  Table  4.1,  i.e.  P  =  49.6  cmVs,  a  =  3.13  cm,  5  =  1 
cm.  The  vortex  generator  was  positioned  at  h/Do  -  1.72  below  the  surface  and  the 
incidence  angle  was  a  =  20*.  The  corresponding  Froude  number  r/(ga3)l/2  o.29. 

Although  double  reconnection  were  observed  at  other  conditions,  these  particular 
conditions  were  selected  because  they  resulted  in  a  highly  reproducible  evolution  of  the 
vortex  rings.  Figure  4.  IS  is  a  cross-section  picture  of  the  vortex  ring  as  it  approaches  the 
surface,  before  any  significant  interaction  with  the  surface  has  occurred. 

4.3.1  Flow  Visualization 

Figure  4.16  (a)  to  (g)  is  a  sequence  of  photographs  showing  the  double 
reconnection  process.  The  photographs  were  obtained  on  different  realizations  of  the  flow 
at  nondimensional  times  Ft/a^  *  4.56,  5.57, 6.58, 7.59,  8.61, 9.62  and  10.63  as 
indicated.  At  each  instant  a  cross-section  view  of  the  underwater  flow  and  a  shadowgraph 
image  of  the  surface  are  presented.  These  were  also  obtained  on  different  realizations  of  the 
flow.  The  two  views  are  arranged  in  the  proper  relative  downstream  position  in  Figure 
4.16.  The  shadowgraph  system  was  imaged  43  cm  above  the  surface  to  obtain  a  better 
quality  image  and  to  increase  the  sensitivity  in  order  to  detea  possible  surface  waves 
generated  during  the  interaction. 

The  vonex  line  reconnection  process  of  the  upper  vortex  is  already  started  at  Ft/a^  = 
4.56  in  Figure  4.16  (a).  This  reconnection  process  is  completed  at  Ft/a^  =  6.58,  Figure 
4.16  (c).  An  interesting  observation  is  the  generation  of  surface  waves  during  the 
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reconnection  process  Ft/a^  =  5.57  to  6.58.  These  waves  which  are  not  present  before  the 
interaction  appear  to  emanate  from  the  interaction  region  and  radiate  away  from  it  (cf. 

Figure  4.16  (b)  and  (c)).  Note  that  there  is  no  dye  associated  with  the  upper  vortex  core  left 
on  the  visualized  plane  after  the  reconnection  process. 

The  vortex  line  reconnection  process  of  the  lower  core  of  the  ring  occurs  between 
Ft/a^  =  7.59  to  9.62  (Figure  4.16  (d)  to  (0).  As  the  vortex  line  reconnection  process 
develops  dye  in  the  cross-section  plane  is  removed.  Also  two  additional  dark  spots  appear 
on  the  surface  signature  indicating  the  presence  of  vortex  lines  terminating  at  the  surface. 
After  the  double  reconnection  process  is  completed.  Figure  4.16  (f)  and  (g),  the  upper  two 
dark  spots  on  the  surface  are  interconnected  by  underwater  vonex  lines.  Similarly  the  lower 
two  dark  spots  on  the  surface  are  also  interconnected  by  vortex  lines  below  the  surface. 

The  two  sets  of  dark  spots  propagate  away  from  each  other  as  shown  in  Figure  (0  and  (g). 
Again  during  the  second  reconnection  process.  Figure  4.16  (f),  there  are  surface  waves 
generated  although  the  propagation  pattern  is  more  complicated  than  in  the  first 
reconnecdon.  In  contrast  to  the  frrst  reconnection,  some  dye  is  left  on  the  visualized  cross- 
section  after  the  second  reconnecdon.  Further  the  cross-section  view  of  Figure  4. 16  (g) 
shows  evidence  of  counterclockwise  vorticity  to  the  left  side  of  the  remaining  dye.  This 
counterclockwise  vorticity  must  be  left  behind  by  the  lower  vortex  core.  This  vorticity  must 
interconnect  the  diver^ng  pair  of  half  v(»tex  rings.  The  total  amount  of  vorticity  in  this 
vortex  appears  to  be  small. 

Detailed  data  on  the  vortex  core  evolution  of  the  double  reconnection  process  were 
measured  from  the  video  recordings.  Simultaneous  visualization  of  the  underwater  flow 
and  surface  deformation  were  obtained  using  the  technique  described  in  Sectioi.  2.3.  Figure 
4. 17  is  a  plot  of  the  downstream  location  of  the  vortex  ring  cores  and  the  distance  of  the 
vc»tex  ring  cores  to  the  free  surface  normalized  by  the  irutial  vortex  ring  diameter,  X/a  and 
Y/a  respectively,  as  a  function  of  the  nondimensional  time  (Ft/a^).  The  speed  of 
propagation  of  the  upper  vortex  cote  is  essentially  constant  until  it  disappears  from  the 
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cross-section.  The  lower  vortex  core  also  moves  with  constant  speed  during  most  of  the 
interaction.  Only  in  the  later  stages  of  the  interaction,  it  slows  down  significantly.  A  least 
square  fit  to  the  data  gives  a  nondimensional  propagation  speed  Uta/r=0.477. 

The  trajectories  of  the  vortex  cores  on  the  free  surface  in  normalized  coordinates  are 
plotted  in  Figure  4.18.  In  this  figure  the  straight  lines  join  the  locations  of  the  upper  and 
lower  vortex  cotes  on  the  surface  at  the  indicated  nondimensional  dme.  The 
nondimensional  propagation  speed  of  the  upper  vortex  core  at  the  initial  state  after 
reconnection  is  0.532  and  the  propagation  angle  with  respect  to  the  centerline  is  54.3*.  For 
the  lower  vortex  core,  a  nondimensional  propagation  speed  of  0.334  and  a  propagation 
angle  of  20.4*  are  measured  at  opening  time.  For  the  steady  state  downstream  of  X/a  »  6, 
the  measured  nondimensional  speed  and  propagation  angle  are  0.280  and  21.2* 
respectively.  After  the  interaction  the  diameter  of  the  resulting  half  vortex  rings  is  increased 
by  19%  compared  to  the  initial  diameter.  The  propagation  speed  is  reduced  by  about  38%. 
This  is  due  to  the  loss  of  vorticity  during  the  interaction  as  well  as  increased  vortex 
diameter. 

4.3.2  Free  Surface  Velocity  Field 

The  velocity  field  at  the  surface  was  measured  using  the  Laser  Speckle  Photography 
technique  discussed  earlier.  As  described  earlier  the  optical  system  used  was  a 
shadowgraph  system  with  the  photographic  plate  positioned  1  cm  above  the  surface  (See 
Figure  3.23).  The  specklegrams  obtained  in  this  manner  contain  not  only  the  panicle 
displacement  infom?ation,  but  also  they  provide  a  low  sensitivity  shadowgraph 
visualization  of  the  surface.  Although  the  simultaneous  visualization  provides  very  useful . 
information  on  the  intoaction  process,  the  intensity  nonuniformity  caused  by  surface 
deformation  prevented  recording  of  the  speckle  pattern  in  several  regions  of  the  flow.  In 
these  experiments  fluid  in  the  generator  was  not  seeded  with  panicles.  Therefore  regions  on 
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the  surface  where  the  fluid  originated  in  the  generator  were  not  seeded  and  consequently  the 
velocity  could  not  be  measured  there. 

Specklegrams  of  the  interaction  were  obtained  at  nondimensional  times  ft/a^  = 

4.56,  5.57, 7.09,  8.10,  9.62  and  10.63.  Contact  prints  of  these  specklegrams  are 
reproduced  in  Rgure  4.19  (a)  to  (f)  respectively.  The  vortex  ring  motion  is  from  right  to 
left.  In  Figure  4.19  (a)  at  Ft/a^  =  4.56  a  curved  dark  region  indicates  a  surface  depression 
caused  by  the  underwater  vortical  flow.  This  region  is  preceded  by  a  wide  region  of  surface 
elevation  and  followed  by  a  much  narrower  elevated  region.  In  Figure  4. 19  (b)  at  Ft/a^  = 
5.57  the  first  reconnection  process  has  been  completed  as  indicated  by  the  presence  of  two 
distinct  daiic  spots  on  either  side  of  the  centerline.  The  dark  spots  are  connected  by  a 
sharply  defined  line.  This  line  separates  the  microballoon  containing  region  of  the  surface 
from  the  clean  region  of  the  surface  produced  by  upwelling  of  fluid  from  below  the 
surface.  A  similar  but  more  diffused  transition  from  clean  surface  to  mictoballoon 
containing  surface  is  observed  on  the  upstream  side  of  the  ring.  The  same  general  features 
are  observed  in  Figure  4. 19  (c)  at  Ff/a^  =  7.09,  the  main  difference  being  the  increased 
distance  between  the  cores.  Figure  4.19  (d)  and  (e)  obtained  for  Ft/a^  =  8.10  and  9.62 
respectively  show  the  surface  deformation  as  the  lower  vortex  core  approaches  the  surface. 
Comparison  with  Figure  4.19  (a)  shows  the  same  local  features  except  that  in  this  case  the 
surface  deformation  region  is  curved  downstream  instead  of  upstream.  Also,  in  these  two 
images  the  transition  from  fluid  with  particles  to  fluid  without  particles  upstream  of  the 
interaction  region  is  sharp  while  downstream  it  is  more  diffused.  Finally  in  Figure  4.19  (0 
at  Ft/a^  s  10.63  when  the  double  reconnection  process  is  completed  two  additional  dark 
spots  are  observed  on  the  surface. 

Vector  plots  of  the  surface  velocity  field  measured  on  the  specklegrams  in  Figure 
4. 19  (a)  to  (f)  are  given  in  Figure  4.20  (a)  to  (0>  As  expiected  the  velocity  could  not  be 
measured  in  very  daric  or  very  bright  regions  of  the  specklegrams  as  well  as  in  regions 
where  there  are  no  particles.  Unfortunately  these  regions  include  the  vonex  cores 
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themselves.  Each  specklegiam  was  interrogated  in  a  square  region  8  x  8  cm  at  0.25  cm 
intervals.  The  velocity  scale  is  the  same  in  all  the  plots  and  is  shown  in  Figure  4.20  (a).  For 
Ft/a^  =  4.56  to  8.10,  Figure  4.20  (a)  to  (d),  the  vector  plots  show  the  velocity  field  induced 
by  two  vortex  cores.  At  Ft/a^  =  9.62  and  10.63  the  measurements  show  the  induced 
velocity  field  of  the  counter  rotadng  vortices.  The  presence  of  this  vonex  results  ia  a 
stagnation  point  and  reverse  flow  in  a  region  on  the  symmetry  plane. 

4.3.3  Deformation  Rate  Field  Analysis 

From  the  measurements  of  velocity  field  at  the  surface  the  components  of  the 
deformation  rate  tensor  in  the  plane  of  the  free  surface  and  the  component  of  the  vorticity 
normal  to  the  h’ee  surface  were  calculated  using  the  methods  described  in  Section  3.3.  In 
addition,  the  component  of  the  deformation  rate  normal  to  the  surface  was  calculated  from 
the  continuity  equation 
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(4.9) 


where  dw/dz  is  gradient  of  the  velocity  component  normal  to  the  surface  in  the  direction 
normal  to  the  surface.  A  positive  value  of  €zz  indicates  that  fluid  below  the  surface  and 

close  to  it  moves  away  from  the  surface.  Similarly  a  negative  value  of  £zz  indicates  motion 

of  the  fluid  below  the  surface  toward  the  surface. 

The  results  (or  the  vorticity  field  showed  very  small  amounts  of  vorticity 
throughout  the  measured  region.  This  is  consistent  with  the  fact  that  the  vonicity  is 
concentrated  in  the  vortex  cores  where  measurements  could  not  be  conducted  due  to  large 
surface  deformation. 

The  results  for  the  deformation  rate  are  presented  in  Figure  4.21.  At  each  point  the 
measured  £xx,  £yy  and  £xy  were  used  to  determine  the  principal  axes  and  corresponding 


43 


deformation  rates.  It  should  also  be  noted  that  at  a  free  surface  in  the  absence  of  capillary 
effects  and  neglecting  the  effect  of  air  above  the  surface,  the  components  of  the  deformation 
rate  tensor  on  the  plane  of  the  free  surface  are  zero.  It  follows  that  the  local  normal  to  the 
free  surface  is  a  principal  axis  of  the  deformation  rate.  Thus  for  small  surface  deflections 
£zz  =  £33  the  deformation  rate  along  the  third  principal  axis.  The  results  are  presented  in 

Figure  4.21  in  the  form  of  a  vector  plot,  labeled  (A)  at  each  flow  condition.  At  each  point 
the  two  principal  axes  are  indicated  by  two  perpendicular  vectors  in  the  directions  of  the 
principal  axes.  The  magnitude  of  each  vector  is  proportional  to  the  magnitude  of  the 
corresponding  deformation  rate  value.  An  arrowhead  on  each  vector  is  used  to  indicate  the 
direction  of  the  deformation  rate.  An  arrowhead  pointing  away  from  the  measurement  point 
indicates  an  extensional  (positive)  deformation  rate.  An  arrowhead  pointing  towards  the 
measurement  point  indicates  a  compressional  (negative)  deformation  rate.  Also  shown  in 
Figure  4.21  are  contour  plots  of  the  £zz  deformation  rates,  labeled  (B)  at  each  flow 
condition.  The  contour  values  plotted  are  Ezz  =  -6,  -4,  -2,  -1,  1, 2, 4,  6  with  the  region  for 
£zz  >  1  (downward  flow)  indicated  by  cross  lines  in  each  plot.  General  examination  of  the 
data  shows  that  the  deformation  rate  £xx,  £yy  and  £xy  arc  very  sensitive  to  the  uncenainty 
of  the  velocity  measurements.  The  uncertainty  of  the  £zz  deformation  rate  is  believed  to  be 
slightly  higher  than  the  x-y  plane  components.  At  Ft/a^  =  4.56,  Figure  4.21  (a),  there  is 

stretching  in  the  direction  perpendicular  to  the  plane  of  symmetry  accompanied  by 
upwelling  (£zz<0)  just  upstream  of  the  interaction  region  near  the  symmetry  plane.  There 
is  a  region  of  positive  £zz  and  consequently  downward  flow  behind  the  interaction.  At 
Ft/a^  SI  5.57,  Figure  4.21  (b),  after  the  first  reconnection  the  deformation  rate  shows  the 
same  basic  features  as  at  Ft/a^  js  4.56.  At  Ft/a^  =  7.09,  Figure  4.21  (c),  the  deformation 
rate  pattern  changes  to  downward  motion  upstream  of  the  interaction.  There  is  also  a  region 
of  upwelling  to  the  side  of  the  interaction  in  front  of  the  vonex  core.  At  Ft/a^  =  8.10, 
Figure  4.21  (d),  the  regions  of  upwelling  and  downward  motion  continue  to  extend 
laterally  following  the  lateral  motion  of  vortex  core.  Upstream  of  the  interaction  a  region  of 
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downward  mouon  is  closely  followed  by  a  region  of  upwelling.  The  magnitude  of  the 
defoimadon  rate  is  small  throughout  the  entire  image.  At  the  beginning  of  the  second 
reconnection,  Tt/a^  =  9.62,  Rgure  4.21  (e),  there  are  regions  of  strong  upwelling  and 
lateral  straining  just  upstream  of  the  interaction.  Finally  at  Ft/a^  =  10.63,  Figure  4.21  (f), 
the  contour  plot  of  8zz  shows  a  larger  region  of  upwelling  upstream  of  the  vortex  cores. 

Also  bands  of  downward  motion  followed  by  upwelling  are  observed.  These  are  believed 
to  be  associated  with  surface  waves.  The  wavelength  1  cm  is  in  agreement  with  the 
wavelength  measured  on  the  photographs  in  Figure  4.16  and  also  they  are  of  the  order  of 
vortex  core  size. 

The  combination  of  flow  visualization  pictures.  Figure  4.17  and  4.19,  surface 
velocity  field.  Figure  4.20,  and  deformation  rate  field  provide  a  detailed  characterization  of 
the  vortex  reconnection  process.  If  we  focus  on  the  reconnection  of  the  upper  vortex  core. 
Figure  4.17  (a),  4.19  (a),  4.20  (a)  and  4.21(b)  suggest  a  slightly  different  dynamics  that 
proposed  by  Schatzle  (1987).  At  Ft/a^  »  4.56,  the  surface  velocity  field  and  surface 
visualization  indicate  that  in  a  stationary  frame  of  reference  there  are  two  stagnation  points 
on  the  surface.  One  is  located  in  front  (to  the  left)  of  the  vortex  core  in  the  light  region  on 
the  shadow  images  (Figure  4.17(a)  and  4.19(a))  the  other  is  located  behind  the  vortex  core 
(to  the  right)  in  the  narrow  bright  region  on  the  shadowgraph  behind  the  kidney  shaped 
dark  surface  depression.  These  two  stagnation  points  are  within  the  region  where  velocity 
measurements  could  not  be  obtained  in  the  velocity  field  plot  in  Figure  4.20  (a),  however, 
they  are  consistent  with  the  overall  flow  pattern.  The  kidney  shaped  dark  region  between 
the  two  stagnation  points  indicates  that  the  surface  velocity  there  must  be  large  and  directed 
upstream  to  be  consistent  with  the  cross-sectional  view.  The  deformation  rate  field  further 
shows  that  in  the  neighborhood  of  the  front  stagnation  point  there  is  upwelling,  a  small 
comnressional  deformation  rate  in  the  flow  direction  and  a  large  extensional  deformation 
rate  in  the  direction  perpendicular  to  the  symmetry  plane.  Therefore,  near  this  stagnation 
point  the  component  of  vorticity  parallel  to  the  surface  and  perpendicular  to  the  symmetry 
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plane  is  enhanced  by  stretching  as  described  by  Schatzle  (1987).  In  contrast  the  flow  near 
the  stagnation  point  behind  the  vortex  core  shows  a  downward  motion  which  will  result  in 
enhancement  of  the  vorticity  component  perpendicular  to  the  surface.  The  velocity  vector 
plot  suggests  compressional  deformation  rates  on  the  free  surface.  Due  to  the  limited 
number  of  point  measurements  of  the  deformation  rate  were  not  obtained  in  this  region. 
Such  a  compressional  deformadon  rate  would  result  in  an  exponential  reduction  of  the 
vorticity  components  parallel  to  the  surface. 

It  is  apparent  that  deformation  rate  near  the  stagnation  point  behind  the  vortex  core 
will  result  in  the  generation  of  vorticity  perpendicular  to  the  free  surface.  As  this  vorticity  is 
generated  the  induced  velocity  field  reduces  the  upstream  velocity  between  the  two 
stagnation  points  and  moves  the  stagnation  points  closer  to  each  other  until  they  merge.  At 
Ft/a^  =  5.57  the  surface  visualization  pictures  and  surface  flow  field  indicate  that  the  flow 
has  reached  this  state. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

Perh^s  the  most  striking  result  of  this  investigation  is  the  detailed  documentation 
of  the  double  reconnection  process  occurring  in  the  interaction  of  a  vortex  ring  with  the  free 
surface.  In  this  type  of  interaction  the  initially  axisymmetric  vortex  ring  splits  into  two 
"half  rings  each  consisting  of  semicircular  vortex  lines  beginning  and  terminating  at  the 
surface.  These  two  "half  rings  propagate  along  diverging  paths  parallel  to  the  surface.  It 
was  found  that  surface  waves  are  generated  during  each  of  the  two  reconnection  processes. 
The  phenomenon  was  observed  over  a  range  of  vortex  ring  Reynolds  number  from  2000  to 
7000  and  incidence  angles  from  10*  to  45*.  The  fu^t  reconnection  process,  that  of  the 
upper  core  was  always  observed  up  to  an  incidence  angle  of  45*.  At  low  Reynolds  number, 
short  slug  length  and  small  incidence  angle  the  "half  vortices  formed  by  double 
reconnection  merged  again  after  a  short  distance.  These  results  are  in  general  agreement 
with  observation  of  colliding  vortex  rings  by  Kambe  and  Takao  (1971),  Fohl  and  Turner 
(1975),  Oshima  and  Asaka  (1977),  Schatzle  (1987)  and  Oshima  and  Izutsu  (1988). 

A  systematic  study  of  the  vortex  reconnection  process  indicates  that  the 
reconnection  time  normalized  by  core  parameters  is  independent  of  Reynolds  number.  The 
reconnection  time  does  depend  on  the  incidence  angle  of  the  vortex  ring.  A  detailed  analysis 
of  the  reconnection  process  using  flow  visualization  and  surface  velocity  field  data  indicates 
that  the  deformation  rate  field  in  the  neighborhood  of  an  stagnation  point  formed  behind  the 
vonex  lines  is  responsible  for  the  rapid  decay  of  the  vorticity  component  parallel  to  the  free 


47 


surface  and  the  rapid  increase  by  stretching  of  the  vorticity  component  perpendicular  to  the 
free  surface,  characteristic  of  the  reconnection  process.  Viscous  diffusion  does  not  play  an 
determining  role  on  the  vortex  reconnection. 

At  zero  incidence  angle,  when  the  vortex  ring  moves  initially  parallel  to  the  free 
surface,  the  distance  from  the  vonex  generator  to  the  interaction  point  and  the  extent  of  the 
interaction  were  found  to  depend  only  on  the  normalized  depth  of  the  vortex  generator. 
They  are  independent  of  the  vortex  ring  Reynolds  number. 

Surface  contamination  by  surface  active  agents  was  found  to  strongly  influence  the 
interaction  of  vonex  rings  with  the  free  surface.  At  an  incidence  angle  of  a  =  0*  with  a 
contaminated  free  surface  the  upper  vonex  core  reconnection  process  was  still  observed. 
The  interaction  of  the  lower  vonex  core  resulted  in  significant  generation  of  opposite  sign 
vorticity  and  rebounding  of  the  core  away  from  the  surface.  At  an  incidence  angle  of  a  = 
90*,  normal  incidence,  the  interaction  with  a  highly  contaminated  surface  resulted  in 
significant  opposite  sign  vonicity  generation  at  the  free  surface.  This  vorticity  rolled  up  into 
secondary  vortices  in  a  manner  similar  to  the  interaction  with  a  solid  surface.  With  a  cleaner 
surface  the  results  showed  significantly  reduced  generation  of  opposite  sign  vorticity  and  a 
delay  on  the  formation  of  the  secondary  vortex.  A  more  systematic  study  of  the  effect  of 
surface  contamination  is  necessary  to  clarify  these  interaction  processes.  In  particular  the 
effect  of  surface  contamination  on  the  double  reconnection  of  vortex  rings  at  a  free  surface 
need  to  be  systematically  investigated. 

A  new  technique  based  on  Laser  Speckle  Photography  was  used  to  measure  the 
velocity  field  at  the  surface.  The  main  limitation  of  the  technique  as  implemented  in  this 
study  was  the  effect  of  surface  deformation  and  seeding  which  limited  the  region  where 
velocity  measurements  could  be  obtained.  The  advantage  in  this  configuration  is  that  in 
addition  to  the  velocity  field  a  simultaneous  visualization  of  the  free  surface  deformation  is 
obtained.  A  new  digital  image  analysis  technique  was  developed  to  determine  the  velocity 
from  the  speckle  photographs.  This  technique  was  tested  and  found  to  have  an  accuracy 
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better  of  1%  in  the  determination  of  the  velocity  field.  Bom  the  velocity  field  the 
components  of  the  vorticity  and  the  deformation  rate  tensor  normal  and  in  the  plane  of  the 
free  surface,  respectively,  were  determined.  The  principal  axes  and  associated  principal 
values  of  the  deformadon  rate  tensor  were  also  determined  as  well  as  the  deformation  rate 
in  the  direcdon  nrxmal  to  the  surface. 


FIGURES 
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Figure  1 . 1  Schematic  diagram  of  the  flow  geometry. 

r,  Circulation;  a.  Ring  diameter,  5,  Core  size; 
(X,  Incidence  angle;  h.  Depth. 
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Figure  2. 1  Schematic  diagram  of  the  facility. 


Figure  2.3  Detail  of  the  vortex  generator  assembly. 


Water 


Vortex  Generator 


Figure  2.4  Vortex  generator  apparatus. 
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Figure  2.5  Block  diagram  of  votex  generator  driving  circuit. 


Figure  2.6  Schematic  diagram  of  simultaneous  flow  visualization. 


Image  Plane  Lens 
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Figure  2.7  Schematic  diagram  of  imaging  system. 


Figure  2.8  Block  diagram  of  hot  film  velocity  measurement  system. 


le  rotating  disk  expenment. 


Figure  3.2  Photograph  of  a  typical  fringe  image. 
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Figure  3.3  Schematic  diagram  of  the  digital  image  processing  system. 


le  Feedback  Processor. 
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Figure  3.5  Power  spectral  density  normalized  to  unit  area  vs.  wavenumber  of  fringe 
image. _ ,  Horizontal  spectrum; . Vertical  spectrum. 


Figure  3.7  Photograph  of  line  averaged  image  for  fringe  image  in  Figure 


Figure  3.9  Photograph  of  image  after  background  subtraction 


wavenumber  ( 1/ lixel) 


Figure  3. 10  Power  spectral  density  vs.  wavenumber  of  the  image  after  background 
subtraction. _ ,  Horizontal  spectrum; . Vertical  spectrum. 


offsct(pixels) 


Figure  3.13  Autocorrelation  vs.  offset  for  image  after  contrast  enhancement. 

_ ,  Horizontal  autocorrelation; . .  Vertical  autocorrelation. 
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1 

I  Figure  3.16  Flow  chart  of  digital  image  processing  program. 


75 


10  m/s 


Figure  3.17  Vector  plot  of  the  velocity  field  measured  in  the  rotadng  disk  experiment. 


minor 


I 


Image  Plane  Lens 


Figure  3. 1 8  Schematic  diagram  of  CVL  free  surface  imaging  system. 
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Figure  3.19  Photograph  of  the  microballoons.  Scale:  full  width  1mm. 


Figure  3.20  Fringe  image  obtained  with  the  surface  imaging  optical  system. 
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Figure  3.2 1  Vector  plot  of  the  velocity  field  t  the  free  surface.  Jet  depth  h/Do  =  1 .0. 

Jet  exit  velocity  25  cm/s.  Framed  area  corresponds  to  2<x/Do<9  and  -2<y/Do<2. 


Figure  3.22  Contour  plot  of  the  vorticity  at  the  free  surface  derived  from  the  results  plotted  in  Figure  3.21. 
Contour  values  plotted  <oDo/Ue  =  -0.8,  -0.5,  -0.3, 0.3, 0.5, 0.8,  1. 
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Figure  3.23  Schematic  diagram  of  CVL  shadowgraph  system  used 
to  obtain  the  speckle  photographs. 
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Figure  4. 1  Velociiy  trace  measured  on  the  centerline  at  X/Do=  5  for  r=  64  cm^/s. 


Position(cm) 


Figure  4.2  Measured  vortex  ring  posidon  as  a  function  of  time  for  several  generator 
operating  conditions. 


Figure  4.3  Measured  Reynolds  number  as  a  function  of  nondimensional  downstream 
position  fcM*  several  generator  operating  conditions. 


Measured  Reynolds  number  after  vortex  ring  formation  vs.  measured  Reynolds  number  at  the  exit. 


■  Re-1200 


(a)  rt/a2  =  4.65 
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Figure  4. 1 1  Sequence  of  photographs  showing  the  interaction  of  a  vortex  ring  with  a  free  surface 
at  an  incidence  angle  a  =  90*.  (A)  Contaminated  surface,  (B)  Cleaner  surface. 
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Figure  4. 1 1  Concluded. 
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Y/Ro 


Figure  4.12  Trajectory  of  vonex  cart  for  the  contaminated  tiee  surface,  (a  =  90*) 

S:  Location  of  primary  vertex  when  secondary  vortex  is  first  observed. 
T:  Location  of  primary  and  secondary  vextex  when  tertiary  vortex  is  first 
observed. 


Figure  4. 14  Nondimensional  reconnection  time  as  a  function  of  Reynolds  number  for 
several  incidence  angles  a  and  depth  h/Do. 
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Figure  4.15  Cross-section  view  photograph  of  a  vortex  ring  before  interaction  with 
a  free  surface.  F/v  =  5000,  o  =  20*. 


1  cm 


(a)  rt/a^  =  4.56 


Figure  4.16  Sequence  of  photographs  showing  the  interaction  of  a  vortex  ring  with 
a  free  surface.  Top  image  :  Shadow  image  of  surface  deformation. 

Bottom  image  :  Cross  section  view.  Dv  =  5000,  a  =  20*. 


Figure  4.16  Continued. 


(c)  Ft/a^  =  6.58 


Figure  4.16  Continued. 


(d)  rt/a2  =  7.59 


Figure  4.16  Continued. 


(e)  ri/a2  =  8.61 


Figure  4. 1 6  Continued. 
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Figure  4.16  Continued. 


Figure  4.16  Concluded. 


Figure  4. 17  Position  of  vortex  cores  in  cross-section  view  as  a  function  of  time, 
r/v*  5000,  a  =  20’. 


Figure  4.18  Trajectory  of  vtxtex  core  on  the  free  surface. 
rAf-5000.a«20*. 


(a)  m2  *4.56 

(Right  edge  Left  edge  X/Do-5.3) 


Figure  4. 19  Sequence  d  speckle  photographs  of  double  reconnection  by  interaction  of 
a  vortex  ring  with  a  free  surface,  TN  »  5(XX),  a  =  20*. 


(b)ri/a2s5.57 

(Right  edge  X/Do*1.3,  Left  edge  X/Do=5.3) 


Figure  4.19  (Continued. 
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(d)  rt/a2  *  8.10 

(Right  edge  ?^o«2.3.  Left  edge  X/Do=6.3) 


Figure  4.19  Continued. 
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■  (c)  rt/a2  =  9.62 

(Right  edge  )^o=2.8.  Left  edge  X/Do=6.8) 
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Figure  4.19  Continued. 
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Figure  4,20  VcctOT  plot  of  a  free  stoface  velocity  field  obtained  from  speckle 
photographs  in  Figure  4.19. 
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Figure  4.20  Continued. 
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Figure  4.20  Continued. 
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Figure  4.20  Continued. 


(0  rt/a2  =  10.63 


Figure  4.20  Concluded. 


Sequence  of  vector  plot  and  contour  plot  of  free  surface  strain  field  obtained  from  velocity  measurements. 
(A)  Vector  plot  of  £i  i  and  £22.  (B)  Contour  plot  of  £33  (Cross  line  region;  £33  >0). 


Figure  4.21  Continued. 
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Figure  4.21  Continued. 


Figure  4.21  Concluded. 
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Table  4. 1  Measurements  of  Vortex  Ring  Characteristics 
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*  r ;  Estimated  from  Figure  4.5. 


APPENDIX  A 


JET  TANK  VORTEX  RING  GENERATOR 


Experiments  were  conducted  using  a  jet  tank  vortex  ring  generator.  These  facility  is 
a  modification  of  the  jet  tank  used  by  Madnia  (1989)  to  generate  vortex  rings  parallel  to  a 
free  surface.  The  free  surface  tank  and  optical  system  arrangement  were  the  same  as 
described  in  Figure  2.6  for  the  simultaneous  flow  visualizations  of  the  underwater  flow  and 
ri-ee  surface  deformation.  The  set-up  for  removing  the  surface  contaminants  was  not 
installed  for  this  experiment  and  the  small  jet  tank  (53x76x30  cm)  instead  of  the  vonex 
generator  system  was  used  to  generate  vonex  rings.  The  small  jet  tank  which  was 
constructed  from  PVC  and  Acrylic  plates  is  shown  in  Figure  A.  1 .  The  jet  exit  orifice, 
which  had  a  circular-arc- shaped  nozzle  with  a  radius  equal  to  the  wall  thickness,  had  a 
circular  cross-section  2,54  cm  in  diamerer.  The  vortex  rings  were  produced  by  a  short 
duration  pressure  pulse  applied  to  the  small  jet  tank  using  the  solenoid  valve  (ASCO 
8262A214)  and  timing  circuit  A  hot  film  prove  at  the  jet  exit  plane  was  used  to  measure  the 
velocity  trace  with  respect  to  the  time  to  determine  the  vonex  ring  circulation.  Because  of 
the  very  large  mass  of  water  in  the  jet  tank  the  velocity  trace  at  the  exit  showed  very  long 
duration  of  water  flow.  This  was  considered  inadequate.  For  more  systematic  studies  the 
vortex  generator  described  in  Section  2.2  was  designed  and  built  to  avoid  this  problem. 
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Figure  A.  1  Schematic  diagram  of  a  small  jet  tank  vonex  ring 
generation  facility. 
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